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NOMENCLATURE 


A  Constant  depending  on  area,  geometry  and  surface  nature  of  slider 
<i,  i  Contact  dimensions 

c  Ratio  of  area  to  weight  jf  slider 
f  Frictional  force 

F.  Force  analogous  to  heat  conducted  into  ice 
F  Force  analogous  to  heat  required  to  melt  surface  of  ice 
F^  Force  analogous  to  heat  conducted  into  slider 
//  Convection  coefficient 
/;  Thickness  of  meltwater  film 
k  Thermal  conductivity 
L  Latent  heat  of  fusion 
N  Load  on  slider 

11  Geometrical  factor 

Q  Rate  of  heat  generation 
q  Heat  flux 

r  Radius  of  contacts 

R  Heat  flow  ratio 

Tm  Melting  temperature 
T()  Ambient  temperature 
t  Time 

v  Velocity 

K  Thermal  diffusivity 
l-i  Coefficient  of  kinetic  friction 
Hm  Energy  dissipated  by  melting 

H  Dynamic  viscosity  of  water 

p  Density 

Subscripts: 

/  Ice  or  snow 

m  Melting 

s  Slider 

w  Water 


IV 


Thermal  Response  of  Downhill  Skis 


GUY C  WARREN ,  SAMUEL C.  COLBECK  AND  FRAr  JC1S  E.  KENNEDY 


INTRODUCTION 


interactive  materials,  and  the  normal  force  N,  or 


The  low  friction  of  snow  and  ice  has  long  been 
recognized  as  one  of  the  distinguishing  character¬ 
istics  of  the  cold  regions.  The  low  friction  is  thought 
to  result  from  the  frictional  heating  produced  by 
the  slider.  The  heat  generated  appears  to  create  a 
thin  meltwater  layer  that  provides  lubrication  and 
reduces  the  area  of  solid-to-solid  contact  between 
the  snow  and  the  gliding  surface  of  the  slider.  The 
amount  of  hea1  generated  at  the  frictional  interface 
accounts  for  more  ilion  9S%  of  c\!l  tilt?  frictional 
energy  losses,  ana  is  very  sensitive  to  the  ambient 
conditions.  The  propagation  of  this  heat  from  the 
frictional  interface  should  depend  on  both  the  am¬ 
bient  temperature  and  the  thermal  characterise  s 
of  the  slider.  Using  both  field  measurements  and  a 
finite  element  model,  we  will  show  that  the  thermal 
characteristics  of  downhill  skis  have  a  significant 
effect  on  the  temperature  at  the  ski  base.  Other 
sliders  should  have  similar  responses. 

Before  we  describe  the  temperature  fields  in 
skis,  we  review  some  of  what  is  already  known 
about  snow  and  ice  friction  and  then  look  at  some 
of  the  most  important  work  done  to  date.  Then  the 
results  of  both  field  measurements  an  '  finite  ele¬ 
ment  simulations  will  be  given. 


Thoughts  on  kinetic  friction 

In  an  effort  to  optimize  sliding  on  snow,  it  is  im¬ 
portant  to  havean  understandingof  thedissipative 
nature  of  fi  iction.Thecoefficient  of  kinetic  friction 
is  defined  as  the  ratio  of  the  frictional  force/ to  the 
normal  force  N  and  is  referred  to  as  Amon ton's 
Law.  When  two  surfaces  rub  together,  the  me¬ 
chanical  energy  loss  is  dissipated  in  a  number  of 
forms,  including  noise  and  heat.  Heat  is  the  dom¬ 
inant  mode  of  dissipation,  so  its  rate  Q  is  given 
approximately  by  the  product  of  the  coefficient  of 
kinetic  friction  p,  the  relative  speed  v  between  the 


Q  =  [ivN  =  vf. 


(1) 


Under  some  conditions,  diis  heat  production 
can  cause  melting  of  one  or  both  of  the  materials  at 
the  frictional  interface;  in  fact,  frictional  melting  is 
of  interest  to;  a  variety  of  materials,  not  just  snow 
and  ice.  When  melting  occurs,  the  resulting  liquid 
acts  as  a  lubricant  that  reduces  the  area  of  solid-to- 
solid  contact  between  the  two  surfaces.  The  reduc- 

in  -V  V  sK.x  tnloefx/'ia  ir  .TxniMa. 
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dent  on  the  thickness  and  viscosity  of  this  liquid 
layer.  According  to  tlu  theory  of  Co’beck  (1988), 
when  an  object  slides  over  snow,  the  sliding  action 
between  the  slider  and  the  snow  can  result  in  melt¬ 
ing  of  the  snow  and  the  int  oduction  of  a  layer  of 
meltwater  that  reduces  the  kinetic  friction  at  the 
interface.  The  meltwater  can  be  sustained  at  the  in¬ 
terface  just  by  the  heat  generated  from  the  shear¬ 
ing  of  the  water.  Thus,  the  sliding  energy  is  used  in 
a  process  that  minimizes  the  friction. 

There  appears  to  be  a  number  of  parameters 
that  affect  the  rate  of  heat  generation  at  the  interface. 
While  the  speed  is  more  or  less  imposed,  the  fric¬ 
tional  force  is  influenced  by  a  number  of  parame¬ 
ters,  including  the  ambient  temperature,  snow 
type  and  the  thermal  conductivity  and  surface 
characteristics  of  the  slider.  Controlled  measure¬ 
ments  of  the  interface  temperature  would  provide 
a  test  of  theextent  to  which  these  parameters  affect 
the  thermal  response  and  therefore  the  sliding  re¬ 
sistance  of  materials  on  snow. 


A  review  of  previous  investigations 

The  idea  that  a  meltwater  layer  is  responsible 
for  the  low  friction  of  snow  and  ice  has  not  always 
been  accepted.  A  study  of  existing  theories  reveals 
that,  although  a  majority  of  the  previous  investi¬ 
gators  were  proponents  of  the  meltwater  theory. 


one  investigator  argues  that  alternative  processes 
other  than  meltvvaterlubrication  might  be  respon¬ 
sible  for  the  lov'  friction  observed.  The  previous 
work  is  worth  looking  at  in  some  detail. 


Bowden  and  Hughes  (1939) 

This  study  laid  the  foundation  for  much  future 
work  on  ski  and  snow  friction.  Bowden  and 
Hughes'  investigation  was  primarily  concerned 
with  the  frictional  characteristics  of  ice,  although 
the  results  apply,  at  least  qualitatively,  to  snow  as 
well  They  showed  that  frictional  heating,  not 
pressure  melting,  is  the  cause  of  melting  of  the 
snow  surface.  This  conclusion  is  easily  demon¬ 
strated  by  calculating  thedepression  of  the  melting 
temperature  from  the  applied  pressure  for  either 
ice  skates  or  skis.  Their  calculation,  however,  was 


made  under  the  assumption  that  t 1  e  entire  base  of 
the  ski  was  in  contact  with  the  snow.  The  real  area 


of  contact  has  been  estimated  by  Kuroiwa  (1977)  as 
3.8%  of  the  apparent  area  of  contact.  If  this  estimate 
is  accurate,  the  ambient  snow  temperature  would 
still  have  to  be  above  -0.004°C  to  allow  pressure 
melting  of  the  snow  at  the  ski/snow  interface. 
Bowden  and  Hughes  also  suggested  that  pressure 
melting  could  not  be  responsible  for  the  low  friction 
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coefficients  differ  greatly. 

The  frictional  heat  is  concentrated  at  the  points 
of  contact  and,  although  much  of  this  heat  is  lost  by 
conduction  into  the  surrounding  ice  and  ski,  a 
sufficient  amount  can  be  retained  to  cause  local 
melting  of  the  snow  crystals  at  the  contact  points. 
This  i  1  lustrated  in  Figure  1 ,  which  shows  a  melt- 
cap  on  a  snow  grain  that  was  removed  by  us  from 
a  snow  surface  after  a  slider  passed  over  it. 

Making  further  rough  calculations,  Bowden  and 


T igurc  2 .  Crystal  from  snow  surface 
with  melt  cap  (T  =  -4°C). 


Hughes  determined  that  sufficient  frictional  heat 
is  generated  duung  skiing  to  raise  the  temperature 
of  the  interface  from  -20  to  0°C  and  that  the  melt¬ 
water  layer  thickness  at  this  ambient  temperature 
is  53  pm  for  a  real  contact  area  of  1 0  mm2.  They  a:  o 
concluded  that  the  extent  of  melting  is  dependent 
on  the  ambient  temperature.  At  very  low 
temperatures,  insufficient  heat  is  generated  tocause 
melting  and  snow  exhibits  the  same  frictional 
characteristics  as  other  unlubricated  materials.  One 
important  point  that  Bowden  and  Hughes  did  not 
discuss  is  that,  although  insufficient  heat  is  gener¬ 
ated  to  initiate  melting  at  these  low  temperatures, 
the  amount  of  heat  generated  is  far  greater  because 
of  the  higher  friction  »hat  exists  in  the  absence  of  a 
thin  meltwater  layer. 

In  their  study  they  directly  measured  the  coef¬ 
ficient  of  kinetic  friction  by  placing  a  smah  s'ider 
on  an  ice-covered  turntable.  When  the  turntable 
was  set  in  motion,  they  were  able  to  measure  the 
frictional  force  by  recording  the  deflection  of  a 
spring.  With  this  method,  they  made  the  following 
qualitative  observations:  Over  a  certain  range  of 
loads,  the  coefficient  of  kinetic  friction  remained 


constant  but  at  extremely  heavy  loads,  the  coeffi¬ 
cient  of  kineuc  friction  decreased.  This  was  a  de- 
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Bowden  and  Hughes,  can  be  attributed  to  the  thick 
meltwater  layer.  They  suggested  that  the  presence 
of  a  thick  meltwater  layer  increased  the  friction, 
due  to  viscous  shear  and  surface  tension,  but  that 
the  coefficient  of  kinetic  friction  was  unaffected  by 
a  change  in  the  apparent  area  of  contact.  Bowden 
and  Hughes  hypothesized  that  speed  influenced 
only  the  temperature  increase  of  the  slider  and  not 
the  temperature  increase  of  the  snow.  They  attrib¬ 
uted  this  to  a  slider  continuously  receiving  heat 
during  sliding,  while  a  section  of  snow  only  re¬ 
ceived  heat  while  the  slider  passed  over  it. 

Bowden  and  Hughes  drew  some  very  interest¬ 
ing  conclusions  abou  t  the  effect  of  the  thermal  con¬ 
ductivity  of  the  slider  on  the  coefficient  of  kinetic 
friction.  Themost  obvious  was  that  a  slider  of  poor 
thermal  conductivity  would  facilitate  more  melting 
because  of  the  greater  retention  of  heat  at  the  fric¬ 
tional  interface.  The  effect  of  thermal  conductivity 
on  melting  was  thought  to  be  more  pronounced  at 
lower  temperatures.  In  an  effort  to  determine  the 
amount  of  melting  occurring  as  a  result  of  frictional 
heating,  they  investigated  the  electrical  conduc¬ 
tivity  between  the  two  counter  surfaces  at  tem¬ 
peratures  close  to  0°C.  They  discovered  that  the 
electrical  conductivity  increased  as  the  slider  was 
set  in  motion,  suggesting  the  formation  of  a 
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continuous  meltwater  layer.  At  lower  tempera¬ 
tures,  tlu>i  igh,  they  observed  smaller  values  of  the 
electrical  conductivity,  suggesting  that  the  melt¬ 
water  layer  was  not  continuous  (i.e.,  melting  oc¬ 
curred  at  localized  points).  Bowden  and  Hughes 
also  observed  the  effect  of  waxes  on  the  bases  of 
skis  and  found  that  the  substitution  of  a  hydro- 
phobic  surface  for  a  hydrophilic  surface  reduced 
the  kinetic  friction.  This  may  be  attributable  toa  re¬ 
duction  in  surface  tension,  but  an  exact  mechanism 
was  not  proposed.  In  any  case,  the  result  reaffirmed 
the  theory  of  the  meltwater  layer. 

A1  though  most  of  their  conclusions  were  based 
on  the  interaction  between  certain  materials  and 
ice,  they  did  a  few  experiments  using  snow  as  the 
underlying  material.  They  concluded  that  the  gen¬ 
eral  behavior  was  very  similar  to  that  observe,  on 
ice,  except  'hat  the  friction  was  always  higher  on 
snow.  This,  they  remarked,  u\.s  ascribable  to  the 
extra  mechanical  work  required  to  deform  the 
snow  surface.  Although  many  assumptions  were 
made  in  their  analysis,  the  results  they  obtained 
were,  at  least  qualitatively,  significant. 

Klein  (1347) 

Klein  studied  the  sliding  resistance  of  aircraft 
skis.  Although  the  skis  used  m  his  investigation 
were  considerably  larger  m  size  and  unit  loading 
(load /total  area),  bis  experimental  observations 
were  almost  identical  to  those  of  Bowden  and 
Hughes.  Klein  explained  his  results  in  terms  of  the 
meltwater  lubrication  theory.  He  suggested  that 
the  layer  of  water  under  the  ski  was  not  ot  uniform 
thickness,  and  that  at  the  front  of  the  ski,  which 
first  came  into  contact  with  the  undisturbed  snow, 
the  surface  was  dry. 

Klein  proposed  that  three  types  of  sliding  resis¬ 
tance  existed  at  the  base  of  a  ski:  solid  friction  at  the 
front  of  theslider,  viscous  resistance  due  to  shearing 
of  the  meltwater  layer  and  surface  tension  effects. 
Each  component  would  depend  on  the  tempera¬ 
ture,  snow  structure,  unit  loading  and  material  in 
contact  with  thesnow.  Klein  assumed  theexistence 
of  solid -to-solid  friction  at  the  front  of  the  ski 
through  evidence  of  wear,  which  ho  noticed  only 
occurred  at  the  front  of  the  ski. 

He  also  made  observations  of  the  real  area  of 
contact.  Previous  values  had  only  been  estimated, 
but  observations  of  water  droplets  in  contact  with 
the  ski  yielded  values  of  20  and  50%  real  contact 
areas  at  unit  loadings  of  1.38xlO(’Pa  (200  lb/ft’) 
and  3.45x1 01’  Pa  (500  lb/ ft-)  respectively.  If  these 
values  are  accurate,  for  a  pressure  of  4.14xl0‘!  Pa 
(60  lb/ft-)  (a  common  value  for  a  skier),  a  6%  real 


contact  area  might  be  assumed .  This  is  comparable 
to  the  estimate  made  by  Kuroiwa  (1 977).  However, 
Klein  did  not  state  under  what  conditions  these 
observations  were  made,  and  the  amount  of  melt¬ 
water  at  the  base  of  the  ski  dependson  the  ambient 
conditions  and  the  thermal  characteristics  of  the 
slider. 

Erick>>on  (1343) 

Ericksson  described  the  effects  of  snow  type  on 
the  coefficient  of  friction  of  skis  and  concluded 
that  the  friction  decreased  ns  the  snow  grain  size 
increased.  The  real  area  of  contact  may  have  been 
sensitive  to  the  grain  size,  particularly  if  the  size  of 
the  asperities  on  the  gliding  surface  of  the  slider 
was  compatible  to  the  dimensions  of  the  snow 
grains.  This  result  supports  the  observation  that 
skis  run  h.ster  on  snow  that  had  undergone  meta¬ 
morphosis  (i.e.,  an  increase  in  grain  size).  This 
grain  size  effect  is  explained  by  Colbeck  (1988)  in 
terms  of  the  dynamics  of  the  water  film. 

Mi  Coimti  (1350) 

McConka  rejected  the  theory  of  meltwater  lu¬ 
brication  and  proposed  that  the  low  friction  on 
snow  was  attributable  to  the  adsorption  of  mono- 
lavers  of  water  vanor  that  nrevent  adhesion,  and 
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thusdecrease  the  resistance  to  starting.  This  theory 
was  proposed  as  a  result  of  a  study  of  the  "diy  s'ip- 
periness"  of  graphite  and  magnesium.  Thes*'  ma¬ 
terials  produce,  underpressure,  plate-iike  crystals 
during  atmospheric  exposure  that  result  in  vapor- 
film  lubrication  during  sliding.  Although  this  the 
ory  indicated  that  another  form  of  lubrication 
might  control  the  basal  friction,  McConica  agreed 
that  ski  geometry  and  composition,  temperature 
and  condition  of  the  snow,  speed  and  loading  af¬ 
fected  the  sliding  resistance. 

McConica  employed  thesameexperimentalpro- 
cedure  as  Bowden  and  Hughes,  but  performed  his 
tests  with  a  greater  degree  of  control  over  the  para¬ 
meters  aforementioned.  Although  the  majority  of 
his  conclusions  were  based  on  theoretical  cal¬ 
culations,  the  experiments  that  he  conducted  did 
yield  some  interesting  results.  He  discovered  that 
a  magnesium  (highly  conductive)  ski  had  a  lower 
resistance  to  sliding  than  a  well-lacquered  wood 
(poorly  conductive)  ski  at  temperatures  well  below 
the  melting  point  of  the  snow.  T urther,  heat  trans¬ 
mission  calculations  he  made  showed  that  the  for¬ 
mation  of  a  meltwater  layer  at  these  temperatures 
was  doubtful.  He  believed  this  suggested  that  the 
retention  of  heat  at  the  base  of  the  ski  did  not  in¬ 
fluence  the  coefficient  of  kinetic  friction  and  that 
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the  low  friction  observed  must  have  been  due  to 
some  other  lubricating  mechanism.  McConicn 
predicted  that  the  base  of  the  ski  could  not  be  at 
any  temperature  other  than  the  ambient 
temperature  of  the  snow  when  the  ski  was  in 
motion.  As  shown  later,  this  idea  is  certainly 
incorrect  and  leads  us  to  question  his  other 
conclusions. 

McConica's  vapor  theory  was  supported  by  the 
fact  that  magnesium  is  a  livdrophilic  material, 
which  would  create  greater  surface  tension  than  a 
hydrophobic  material  (like  lacquered  wood)  in  the 
presence  of  a  meltwater  layer.  In  addition,  a 
comparison  between  magnesium  and  aluminum 
showed  that  the  magnesium  ski  had  a  lower 
coefficient  of  kinetic  friction  under  all  conditions 
because  the  film  that  forms  on  aluminum  is 
amorphous  and  is  not  susceptible  to  vapor  film 
lubrication.  These  results  should  be  tested  further, 
especially  in  view  of  magnesium's  chemical 
reactivity  with  water.  Additional  observations  by 
McConica  suggested  that  kinetic  friction  decreased 
as  load  increased,  and  kinetic  friction  decreased  as 
speed  increased. 

Hou’ticu  (1953) 

Bowden  conducted  tests  with  real  snow  skis  to 
determine  theeffect  of  ski  base  finish,  temperature, 
speed  and  snow  hardness  on  friction.  He 
determined  that  as  the  speed  increased,  the 
coefficient  of  friction  was  reduced  because  of  the 
increase  in  melting.  Also,  as  the  temperature 
decreased,  the  coeffich"  increased  because  of  the 
lack  of  melting.  According  to  Bowden,  a  coating  of 
Teflon  on  the  base  ot  the  slider  resulted  in  lower 
friction  than  any  waxed  base.  This,  he  concluded, 
was  ascribable  to  the  high  hydrophobicity  of  Teflon. 
In  tests  using  waxes  as  a  base  material,  he  found 
that  the  snow  was  capable  of  penetrating  the  wax 
layer,  resulting  in  a  higher  coefficient  of  friction 
than  for  the  harder  Teflon. 

Bowden's  study  further  supported  his  theory 
that  localized  surface  melting,  caused  by  frictional 
heating,  was  responsible  for  the  low  friction  of 
sliders  on  snow.  His  analysis  of  the  hydrophobicity 
of  certain  materials  also  explained  the  extent  to 
which  surface  tension  affected  the  ability  of  a  ski  to 
move  forward  in  "wet"  conditions. 

Tufimn  ami  Yosida  (1969) 

The  experimental  study  carried  out  by  Tusima 
and  Yosida  best  replicated  real  skiing  conditions. 
Most  previous  experiments  slid  a  ski  on  a  moving 
turntable,  which  allowed  the  ski  to  pass  over  the 


same  section  of  snow  for  an  indefinite  per  iod .  This 
accumulation  of  heat  in  both  the  ski  and  the  snow 
would  certainly  increase  the  likelihood  of 
meltwater  formation.  In  the  case  of  real  skiing,  the 
snow  is  always  renewed,  which  limits  the  extent  to 
which  the  snow  surface  could  be  heated. 

Applying  this  experimental  technique,  Tusima 
and  Yosida  discovered  that  melting  began  as  soon 
as  forward  motion  began,  but  that  the  rate  of 
melting  reached  a  constant  value  with  time.  They 
theorized  that  the  points  of  contact  were  at  0°C 
after  only  0.1  seconds  of  forward  motion  and  from 
their  experiments  determined  that  the  water  film 
thickness  was  between  1 0  and  50  gnv  This  range  of 
values  of  the  thickness  was  lower  than  those  esti¬ 
mated  by  Bowden  and  Hughes,  possibly  because 
of  the  different  experimental  procedures  used. 

Earns,  Ni/c  and  Chrcsaman  (1976) 

Although  Evans  et  al.  were  primarily  interested 
in  the  kinetic  friction  of  ice,  they  were  the  first 
investigators  to  use  thermocouples  in  their  analysis 
of  temperature  rises  associated  with  frictional 
heating  involving  ice  or  snow.  The  measurements 
they  made  enabled  them  to  predict  the  effect  of 
slider  thermal  conductivity  on  friction.  They  found 
that  fora  copper  slider,  between  40  and  60%  of  the 
heat  generated  at  the  trictional  interface  was 
conducted  away  through  the  copper  and  that  this 
fraction  was  apparently  independent  of  the  air 
temperature  between  -2  and  -15°C.  In  addition, 
the  heat  conduction  through  the  slider  depended 
on  the  speed. 

In  determining  the  coefficient  of  friction,  they 
made  an  analogy  between  the  heat  generated  and 
the  frictional  force.  They  defined  the  total  frictional 
force  at  the  contact  points  as 

f~F+  F.+  F  (2) 

where  Fs  =  force  analogous  to  the  heat  conducted 
into  the  slider 

F  =  force  analogous  to  the  heat  diffused 
into  the  ice 

F  =  force  analogous  to  the  heat  required  to 
melt  the  surface  of  the  ice. 

According  to  Evans  et  al. 

r  _  A  k  s(  7  „,-T0)  (3) 

V 

where  A  =  constant  that  is  dependent  on  the 
contact  area,  geometry’  and  surface 
nature  of  the  slider 
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A  =  thermal  conductivity  of  the  slider 
T  -  melting  temperature 
T(1  =  ambient  temperature  of  the  ice 
v  =  relative  speed  between  the  slider  and 
the  ice. 


Also,  they  derived  the  force  analogous  to  the  heat 
diffused  in  to  the  ice  as 


/'  =  2k{T  -T,.)b 

1  l  m  o  1 


(4) 


where  k  =  thermal  conductivity  of  the  ice 
a  =  length  of  contact  points 
b  -  width  of  contact  points 
k  -  thermal  diffusivity. 


I’utting  these  two  expressions  together  with  eq 
1,  the  coefficient  of  friction  becomes 


Ak  AT.. 


vN 


-Tid+1.74f,(7'1 


N(TKl-)0-5 


+  Mi, 


(5) 


where  ji  represents  die  energy  dissipated  by 
melting. 

According  to  Evans  et  al.,  although  jr  increased 

twifh  1'  if  line  n  ciinnll  rnii  frikt  if  i#\ti  fhu 
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friction  and  could  often  be  neglected.  From  eq  5, 
we  see  that  the  coefficient  of  friction  depends  on 
the  thermal  conductivity  of  the  slider  and  the  ice, 
the  area  and  geometry  of  the  contact  points,  the 
ambient  and  melting  temperature  of  the  ice,  the 
diffusivity  of  the  ice,  the  load  applied,  and  the 
speed  ot  the  slider.  Comparing  the  values  of  friction 
obtained  in  their  experiments  with  the  values  cal¬ 
culated  by  eq  5,  they  found  that  the  values  agreed 
within  a  factor  of  two.  The  discrepancy,  they  con¬ 
cluded,  was  attributable  to  an  overestimation  of 
the  real  contact  area. 


Aml'iU'li  and  Mayr  (19S1 ) 

Ambacli  and  Mayr  were  able  to  determine  the 
thickness  of  the  meltwater  layer  by  measuring  the 
capacitance  at  the  slider/snow  interface.  Their 
measurements  were  made  at  speeds  approaching 
60  km /hr  (approximately  1 7  m/s)  and  they  found 
that  the  thickness  of  the  water  layer  varied  with 
speed,  snow  type,  snow  and  air  temperature,  and 
the  gliding  ski  surface.  They  concluded  that  the 
meltwater  layer  in  skiing  was  5  to  10  pm  thick  and 
that  its  value  increased  with  an  increase  in  both  the 
snow  and  air  temperature,  and  with  an  increase  in 
speed.  In  determining  theeffectofsnow  type,  they 
used  packed  snow  and  deep,  unpacked  snow  and 


found  that  the  malt  water  layer  thickness  was  great¬ 
er  in  the  deep  snow.  They  suggested  that  the  rea¬ 
son  for  this  was  that  on  normal  ski  terrain,  the 
snow  surface  was  bumpy  so  heat  generation  was 
intermittent.  Although  Ambacli  and  Mayr  did  not 
know  the  specific  characteristics  of  the  waxes  they 
applied  to  the  base  of  their  test  skis,  they  discovered 
that  the  meltwater  thickness  was  larger  for  waxes 
whose  temperature  ranges  included  the  ambient 
temperature  of  their  tests.  This  indicated  that  an 
appropriate  wax  appeared  to  increase  the  thickness 
of  the  meltwater  layer,  resulting  in  a  lower  coeffi¬ 
cient  of  kinetic  friction. 

Oksanen  (1 983) 

Oksanen  studied  the  friction  and  adhesion 
strength  of  ice.  They  measured  the  coefficients  of 
both  static  ana  kinetic  friction  between  ice  and 
several  materials  as  a  function  of  load,  speed  ar  d 
temperature.  The  apparatus  used  in  this  study  we  s 
a  modified  version  of  a  turntable  that  Keinonen 
had  used  in  his  study  of  friction  between  skis  and 
snow. 

Oksanen  predicted  that  at  low  temperatures 
(less  than  -10°C),  tin*  coefficient  of  kinetic  friction 
would  vary  as  ;  11 ,  where  v  is  the  relative  speed 

between  the  material  and  the  ire.  This  ,  he  theorized, 
was  ascribable  to  the  •self-balancing  ot  the  thin 
water  layer.  According  to  Oksanen,  if  the  thickness 
of  the  meltwater  layer  decreased,  the  increased 
frictional  heating  would  create  more  water,  and  if 
the  meltwater  layer  thickness  increased,  the  re¬ 
duction  in  frictional  heating  would  cause  a  drop  in 
temperature  at  the  contact  below  the  melting  point 
of  the  water.  So  at  equilibrium,  the  temperature  at 
the  contact  surface  was  at  the  melting  point  of 
water  and  the  heat  produced  by  the  friction  was 
equal  to  the  heat  conducted  into  the  two  solids.  At 
higher  temperatures  (close  to  0°C),  Oksanen  deter¬ 
mined  that  his  equations  reduced  to  a  form  that 
showed  that  the  friction  was  determined  by  the 
viscous  shear  and  increased  as  v  '  *.  His  experi¬ 
mental  results  concerning  the  coefficient  of  kinetic 
friction  supported  Iris  theory  to  a  good  approxi¬ 
mation.  The  discrepancies,  he  concluded,  were 
due  to  the  assumption  that  viscous  shearing  was 
the  only  component  of  friction  when  a  meltwater 
layer  was  present. 

Akkek.  Ettlcs  mid  Calabrese  C19S7) 

Akkok  et  al.  developed  a  theory  that  suggested 
that  the  heat  generated  at  the  frictional  interface  at 
-20°C  was  sufficient  to  raise  the  temperature  of  the 
interface  by  40°C.  They  contested  that  this  rise  in 


temperature  was  not  realistic  because  any  water 
that  formed  at  the  contact  points  was  immediately 
removed  by  the  motion  ot  the  slider.  In  light  of  this, 
they  suggested  that  the  interface  could  not  attain  a 
value  any  higher  than  0°C  and  that  the  friction 
must  be  dependent  on  the  ambient  snow  temper¬ 
ature.  From  their  analytical  model,  they  also  sug¬ 
gested  that  the  coefficient  of  kinetic  friction  depend¬ 
ed  on  the  load,  speed  and  thermal  characteristics 
of  the  slider.  In  addition,  Akkok  et  al.  were  able  to 
predict  the  thickness  '  the  meltwater  layer  and 
obtained  values  that  uere  comparatively  small. 
They  predicted  a  thickness  of  only  a  few  rnicro- 
..  .eters,  and  said  that  the  surface  roughness  of  the 
slider  prevented  any  hydrodynamic  effects. 

Colbeck  (1988) 

Colbeck  developed  a  theory  that  predicts  the 
heat  generation  by  an  object  sliding  on  snow,  the 
thickness  of  the  resulting  meltwater  layer,  and  the 
effects  of  temperature,  slider  thermal  conductivity, 
speed  and  snow  type  on  the  coefficient  of  kinetic 
friction.  This  theory  was  developed  to  describe  the 
intermediate  values  of  the  meltwater  layer 
thickness  where  the  viscous  effects  of  the  meltwater 
film  dominate  frictional  forces  arising  from  solid- 
to-soiid  ploughing  or  surface :  ension  forces.  From 
the  mechanisms  controlling  the  lubricated  compo¬ 
nent  of  friction  gwat  the  melting  temperature,  Col¬ 
beck  found  that 


where  c  is  the  ratio  of  area  to  the  weight  of  the 
slider,  and  h  is  the  th'rkness  of  the  water  film.  In 
addition,  Colbeck  calculated  the  square  of  the 
thickness  of  the  water  film  to  be 

fi2=^  (7) 

Pv/ 

where  r  =  radius  of  the  water  film 
p  =  density  of  water  al  0°C 
L  =  latent  heat  of  fusion. 

According  to  Colbeck,  this  implied  that  even  in 
fairly  icy  conditions,  the  thickness  of  the  meltwater 
layer  never  exceeded  about  2  (im  in  snow  skiing. 
The  discrepancy  between  this  valueand  the  values 
obtained  by  previous  investigators  may  have  been 
caut  ed  by  the  different  effects  of  hydrophobic  and 


hydrophilic  surfaces.  His  analysis  further  suggest¬ 
ed  that  r  increased  with  repeated  passes.  If  this  is 
true,  then  from  eq  6  and  7,  the  coefficient  of  lubri¬ 
cating  friction  is  reduced  as  r  increases,  as  Ericksson 
(1949)  observed. 

Colbeck  also  predicted  values  of  the  kinetic  fric¬ 
tion  when  the  ambient  temperature  was  below  the 
melting  temperature  of  the  snow.  It  appeared  that 
the  thermal  characteristics  of  the  slider  and  snow 
greatly  affected  the  frict  ion.  At  lower  temperatures, 
the  dry  friction  coi  p..nent  was  larger  than  the 
lubricated  component  because  of  insufficient  melt¬ 
ing.  Calculations  of  the  heat  generated  by  die  dry 
friction  suggested  that  the  front  of  the  slider  was 
d  ry  over  a  certain  length.  At  the  end  of  this  dry  por¬ 
tion,  the  water  layer  increased  ir  th'ckness  towards 
an  equilibrium  va’ue .  The  length  of  the  d  ry  section 
was  dc-j  mdent  on  the  temperature  and  the  slider 
material. 

Colbeck  also  concluded  that  the  lubricated  fric¬ 
tion  increased  with  i’1  /2at  the  melting  temperature. 
However,  at  lower  temperatures,  this  component 
of  friction  decreased  before  increasing  with  n1  /2  at 
higher  speeds.  In  addition,  Colbeck  concluded 
that  the  thickness  of  the  meltwater  layer  was  highly 
dependent  on  the  size  of  the  contacts,  speed,  tem¬ 
perature  and  thermal  vliaiaeleusliss  of  the  slider, 
but  was  independent  of  the  load  applied.  Unfortu¬ 
nately,  many  assumptions  were  necessary  in  these 
predictions. 

Since  much  of  our  knowledge  of  snow  and  ice 
friction  is  based  on  theory,  it  is  important  to  make 
measurements  that  provide  direct  information 
about  the  processes  responsible  for  the  low  friction. 
While  indirect  measurements  of  the  meltwater 
have  been  made,  direct  measurements  of  the  tem¬ 
perature  field  have  not.  Thus,  we  measured  tem¬ 
peratures  to  examine  the  heat  generation  and  heat 
flow  characteristics  of  downhill  skis.  We  also  gen¬ 
erated  solutions  for  the  temperature  field  using 
the  finite  element  technique  to  show  that  the  heat 
flow  characteristics  of  downhill  skis  can  be  simu¬ 
lated.  This  is  necessary  to  show  how  the  construc¬ 
tion  of  skis  affects  their  heat  flow  characteristics, 
and,  therefore,  their  frictional  behavior. 


FIELD  STUDY 

Measurcrut.iis  uf  rki  temperatures 

To  test  some  c-f  the  ideas  given  above  it  is  very 
desirable  to  measure  the  temperature  of  the  actual 
contact  points  between  a  slider  and  the  snow,  but 
limitations  on  thermocouple  size  and  response 
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times  restrict  us  to  measurements  of  the  bulk  tem¬ 
perature  of  the  slider.  However,  if  the  contact 
points  are  at  the  melting  temperature,  as  has  been 
suggested  many  times,  heat  flow  patterns  up 
through  a  ski  may  indicate  the  extent  of  melting  at 
the  frictional  interface. 

From  the  previous  investigations,  the  heat  gen¬ 
erated  by  friction  appears  to  depend  on  the  speed 
and  thermal  characteristics  of  the  slider,  the  load 
applied  to  the  slider,  the  ambient  temperature  of 
the  snow  and  the  snow  type.  Although  our  ability 
to  control  or  even  measure  these  parameters  was 
rather  limited  in  th  ^  tests  described  below,  we  felt 
that  more  information  could  be  gained  from  the 
higher  speeds  and  natural  snow  conditions  found 
in  skiing  rather  than  on  a  laboratory  test  track.  The 
highest  speeds  attained  in  these  tests  were  about 
18  m/s. 

The  objective  is  to  measure  the  temperature 
rises  in  skis  attributable  to  frictional  heating.  If  the 
amount  of  heat  generated  is  measurable,  it  might 
be  possible  to  analyze  the  effects  of  various  para¬ 
meters  on  kinetic  friction.  These  parameters  in¬ 
clude  speed,  load,  snow  type,  ambient  temperature, 
and  slider  geometry  and  thermal  conductivity.  In 
addition,  measurements  of  the  temperature 
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the  materials  that  would  be  best  at  decreasing  fric- 
cion 

In  developing  an  approach  to  these  measure¬ 
ments,  we  decided  that  thermocouples  installed  in 
a  downhill  ski  could  be  used  to  calculate  the  fric¬ 
tional  heat  generated  at  the  base  and  the  ensuing 
heat  propagation  in  the  ski.  The  thermocouples 
implanted  were  small  so  that  significant  changes 
in  the  structure  of  the  ski  could  be  avoided  and  so 
that  the  time  response  of  the  thermocouples  would 
be  short.  Two  sids  were  selected  for  instrumenta¬ 
tion.  a  Blitz  child's  skiand  a  Rossignol  DH  downhill 
racing  ski.  In  the  Blitz  ski,  the  thermocouples  were 
implanted  by  boring  a  hole  just  wide  enough  for 
them  to  slide  into  their  specified  locations.  Then 
the  holes  were  back-filled  with  an  epoxy  resin  to 
secure  the  thermocouples  in  place,  strengthen  the 
ski  and  minimize  erroneous  heat  flow  effects.  Tire 
DH  ski  was  built  at  the  factory,  where  the  thermo¬ 
couples  were  installed  during  construction. 

Description  of  thermocouples 
and  data  acquisition  system 

It  was  essential  to  employ  a  data  acquisition 
system  that  was  both  durable  and  portable  to 
measure  the  interfacial  temperature  increases  dur¬ 
ing  skiing .  In  view  of  this,  one  of  the  small,  computer 
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Figure  2.  Wiring  configuration  between  data  logger 
and  multiplexer. 


con  trolled  da  ta  loggers  was  chosen .  The  CR 1 0  was 
compact  enough  to  be  carried  by  a  skier  in  a  back¬ 
pack  and  rugged  enough  to  survive  the  physical 
abuse  experienced  in  downhill  skiing.  This  system 
consists  of  an  analog-to-digital  converter  that  re¬ 
ceives  voltage  differentials  from  the  thermocouples 
through  the  wiring  panel  or  the  AM  32  multiplexer 
board .  Th  i  multiplexer  board  was  capable  of  sam¬ 
pling  32  thermocouples  and  was  used  when  more 
than  five  thermocouples  were  to  be  sampled.  The 
wiring  between  the  CR  10  and  the  multiplexer  is 
shown  in  Figure  2. 

The  programs  for  controlling  the  Ck  10  are 
given  in  Appendix  A.  Certain  keystrokes  enabled 
the  skier  to  look  at  recorded  and  real-time  temper¬ 
atures.  The  real-time  temperature  readout  was 
convenient  for  determining  when  the  stationary 
ski  was  in  thermal  equilibrium.  All  of  the  data 
were  stored  in  a  solid-state  storage  module  that 
was  small  enough  to  be  put  into  the  backpack  too. 
The  storage  module  was  capable  of  storing  88,000 
double  accuracy  data  points,  more  than  enough 
storage  space,  even  on  a  day  when  30  or  more  runs 
were  taken.  Including  the  two  6-V  batteries,  the 
backpack  and  its  contents  weighed  about  10  kg. 

Type  T  (copper-constantan)  thermocouples 
were  chosen  because  of  their  5-ps  response  time 
and  small  size.  Although  the  0.5-mm-diameter 
thermocouple  lead  wires  had  advantages,  they 
were  too  fragile  for  this  application.  During  the 
course  of  the  field  tests,  a  n  umber  of  thermocouples 
either  gave  erroneous  values  or  quit  completely. 
Some  of  these  were  replaced  during  the  tests. 

Blitz  ski 

This  ski  was  selected  because  its  homogeneous 
structure  greatly  simplified  thermocouple  em- 
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placement  and  interpretation  of  the  heat  flow  field. 
The  Blitz  ski  was  1.1  in  long,  with  an  approximate 
apparent  contact  area  of  0.05  m~.  Five  thermo¬ 
couples  were  installed  as  illustrated  in  Appendix 
B,  yielding  both  transverse  and  vertical  tempera¬ 
ture  distributions  in  the  ski.  As  shown  in  Appendix 
B,  the  thermocouple  arrangement  was  changed 
after  some  use  of  the  Blitz  ski. 

With  only  five  thermocouples,  the  thermocouple 
lead  wires  could  be  connected  directly  to  the  CR 10 
wiring  panel  located  on  top  of  the  CR  10  data  ac¬ 
quisition  system.  This  increased  the  range  of  sam¬ 
pling  frequencies,  although  for  all  runs  with  the 
Blitz  ski,  a  1-Hz  sampling  frequency  was  used. 

PH  ski 

This  ski  is  used  in  World  Cup  downhill  races 
where  ski  glide  is  very  important.  The  ski  was  2.25 
m  long,  with  a  contact  length  of  approximately  2.1 
m  and  an  apparent  contact  area  of  approximately 
0.16  m".  The  DH  ski  was  constructed  and  \<stru- 
mented  by  the  Racing  Department,  Rossignol  S.A., 
in  France.  Appendix  B  (Fig.  B2)  shows  where  32 
thermocouples  were  placed  in  the  DH  ski  at  loca¬ 
tions  specified  by  us.  The  thermocouples  were 
strategically  located  so  that  longitudinal,  transverse 
and  vertical  temperature  profiles  could  be  record¬ 
ed.  1  his  configuration  yielded  a  three-dimensional 
matrix  of  temperatures  in  the  ski.  In  addition,  ther¬ 
mocouples  were  placed  on  one  of  the  steel  edges. 
To  avoid  erroneous  readings,  the  thermocouples 
on  the  edges  had  to  be  insulated  because  of  the 
electrical  conductivity  of  steel.  In  the  four-place 
designation  of  the  thermocouples,  the  number  in 
the  first  pi  ce  indicates  the  longitudinal  position 
in  the  ski  and  the  letter  in  the  second  place  means 
a  vertical  array  (a),  a  basal  array  (b)  or  an  isolated 
edge  or  basal  thermocouple  (c).  Thelast  two  places 
are  for  position  in  the  vertical  array  (02  to  42),  the 
position  along  the  base  (01  to  05),  or  an  edge  ther¬ 
mocouple  (10)  or  an  isolated  basal  thermocouple 
(02). 

Testing 

Various  downhill  ski  areas  were  used  as  test 
sites.  For  each  test,  air  and  snow  temperature, 
weather  conditions  and  snow  conditions  were  re¬ 
corded  at  the  beginning  and  end  of  each  run.  In 
particular,  the  air  and  snow  temperatures  were 
measured  at  the  top  and  bottom  of  each  run  to  ac¬ 
count  for  the  temperature  variation  with  altitude. 
The  snow  temperature  was  measured  bv  placing 
theski  in  contact  with  the  snow  surface  and  reading 
the  temperature  of  a  basal  thermocouple  once  it 
had  stabilized .  The  air  temperature  was  determined 


by  placing  the  ski  about  0.1m  above  the  snow  sur¬ 
face  and  reading  the  steady-state  temperature  of 
the  same  basal  thermocouple.  Because  of  the  fast 
response  time  of  the  thermocouples  and  their  prox¬ 
imity  to  the  base  of  the  ski,  we  observed  a  steady- 
state  condition  within  about  20  seconds  when  de¬ 
termining  the  snow  and  air  temperatures.  Once 
the  test  ski  and  thermocouple  wires  had  been  se¬ 
cured  to  the  skier,  the  data  acquisition  system  was 
programmed  and  the  ski  was  allowed  to  rest  flat 
on  the  snow  until  a  steady-state  temperature  con¬ 
dition  existed  throughout  the  ski. 

Each  test  consisted  of  a  simple,  but  careful,  ski 
run.  Ideally,  more  control  over  the  variables  would 
have  been  possible;  most  of  the  test  conditions 
were  imposed  on  us  by  the  prevailing  conditions 
and  the  mechanics  of  skiing.  Speed  and  load  are 
highly  variable  because  of  changes  in  microtopog¬ 
raphy,  temperature  is  affected  by  altitude  and 
solar  input,  and  snow  type  is  affected  by  temper¬ 
ature,  recent  snow  falls  and  skier  usage.  We  decided 
that,  as  long  as  the  fluctuations  in  these  parameters 
were  recorded  or  reasonably  consistent  during 
each  run,  the  results  would  be  useful. 

Parameters  affecting  frictional  heating 

Ambient  snow  temperature.  Previous  investigat¬ 
ors  have  stated  that  the  amount  of  factional  heat 
generated  is  dependent  on  the  ambient  snow  tc  m- 
perature.  Runs  were  made  in  which  all  of  the  para¬ 
meters  were  consistent  except  the  ambient  snow 
temperature.  There  was  some  difficulty  in  deter¬ 
mining  this  effect  because,  on  any  given  day,  the 
snow  temperature  did  not  vary  more  than  a  4  or 
5°C  between  the  first  and  last  runs.  However,  this 
small  ambient  temperature  differential  appeared 
to  be  sufficient  to  prove  this  effect,  as  is  shown 
later.  Since  the  upper  limit  of  the  basal  temperature 
should  be  the  melting  temperature,  tests  were  con¬ 
ducted  at  temperatures  just  below  0°C  to  see  what 
temperature  the  base  could  reach. 

Load.  We  tested  the  fluctuation  of  load  at  con¬ 
stant  speed  by  cyclically  weighting  and  unweight¬ 
ing  the  ski  approximately  every  20  seconds,  while 
the  skier  was  being  dragged  up  a  Porna  lift.  Once 
the  Poma  had  attained  constant  speed,  the  skier 
applied  approximately  half  of  his  weight  (N/ 2)  to 
the  lest  ski  and,  10  seconds  later,  applied  all  of  his 
weight  to  the  ski  before  returning  to  half  of  his 
weight  10  seconds  after  that.  This  loading  pattern 
continued  for  the  duration  of  the  uphill  run,  at 
which  point  the  ski  was  stopped  so  that  it  could  at¬ 
tain  a  steady-state  condition  again.  Tire  speed  of 
the  slider  in  the  snow  track  was  about  2  m/s. 

Tr  sts  of  constant  load  during  descent  were  done 
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to  have  all  three  contributing  factors  to  the  rate  o( 
heat  generation  (load,  speed  and  coefficient  of 
friction)  kept  constant.  Since  neither  the  speed  nor 
the  coefficient  of  friction  could  be  measured  dur¬ 
ing  the  runs,  very  simple  te^.s  were  conducted  in 
which  only  the  load  varied  between  consecutive 
runs.  Tests  were  conducted  with  loads  of  approxi¬ 
mately  N/2, 3N/4  and  N  applied  to  the  test  ski. 

Tests  of  turning  at  relatively  low  speeds  were 
conducted  with  the  skier  continuously  turning 
either  left  or  right.  The  cyclic  loading  that  occu  rred 
during  turning  should  show  how  the  rate  of  heat 
generation  is  affected  by  this  loading  pattern. 

Tests  utilizing  an  even  weigh*  distribution  with 
the  skier  skiing  in  the  parallel  position  in  the  tuck 
(or  crouch)  were  done.  This  showed  if  the  skier's 
weight  is  evenly  distributed  across  the  width  of 
the  ski. 

Speed.  Tests  were  conducted  to  determine  the 
temperature  increases  that  result  from  a  difference 
in  speed  between  downhill  and  uphill  runs.  The 
skis  w'ere  kept  parallel  on  a  downhill  run  and  the 
speed  increased  at  an  approximately  constant  ac¬ 
celeration.  This  test  was  repeated  for  a  number  of 
different  loads  ranging  from  .V/2  to  N,  and  on  runs 
of  varying  gradient. 

Snoiu  type.  Tests  were  conducted  on  separate 
da  vs  when  the  snow  tvne  had  dramaticallv 

J  J  i  -f 

changed  because  of  changes  in  precipitation  and 
temperature.  Thus,  tests  were  performed  in  snow 
condibons  that  varied  from  sof*,  fresh,  wet  snow 
to  dry,  hard-packed  snow. 

Thermal  conductivity.  One  of  the  most  important 
parameters  affecting  the  temperature  at  the  inter¬ 
face  is  the  thermal  conductivity  of  the  slider.  This 
effect  was  important  to  investigate  because  it  is 
one  of  the  few  things  n  the  design  of  any  slider  that 
can  be  directly  controlled.  It  was  investigated  by 
measuring  the  temperatures  in  the  vertical  direction 
in  both  skis.  Although  identical  ambient  conditions 
never  occurred  during  testing  of  both  skis  (i.e.,  the 
skis  were  never  used  on  the  same  test  day,  but 
were  taken  out  for  testing  on  different  days  and 
therefore  experienced  different  snow  types  and 
ambient  temperatures),  differences  in  the  vertical 
profiles  could  illustrate  the  extent  to  which  thermal 
conductivity  affects  the  flow  of  heat  from  the  fric¬ 
tional  interface. 

Data  processing 

At  the  end  of  a  test  session,  the  data  were  down  ■ 
loaded  from  the  storage  module  to  the  hard  disk  of 
a  personal  computer  The  data  could  then  be  con¬ 
verted  into  temperature  versus  time  plots.  Before 


graphing,  we  smoothed  the  data  from  both  skis 
with  a  0.5-Hz  filter  to  reduce  the  noise  in  the  signal. 

Results  and  discussion 

Because  of  the  highly  dynamic  conditions  ex¬ 
perienced  in  skiing,  it  was  very’  difficult  to  isolate 
one  parameter  and  determine  its  influence  on  the 
frictic  \1  heat  generation.  This  problem  was  en¬ 
countered  throughout  the  experin  ents  and  can  be 
attributed  to  the  frequent  variations  in  direction, 
load  and  speed  that  were  experienced  even  on  ski 
slopes  that  appeared  smooth  and  of  consistent 
gradient  and  snow  type.  In  addition,  it  was  difficult 
to  quantify  any  of  these  parameters.  Nevertheless, 
some  useful  information  was  produced,  although 
the  conclusions  drawn  from  the  study  are  mostly 
qualitative.  The  temperature  versus  time  plots 
from  the  field  study  shown  in  Appendix  C  are  only 
examples  of  the  data  taken;  many  sinvlar  plots  ex¬ 
ist  but  are  not  shown  because  they  are  redundant. 

Effects  of  ambient  temperature 

There  is  a  definite  trend,  as  shown  in  Figures  Cl 
through  C4,  that  suggests  that  the  base  of  the  test 
ski  reaches  a  state  of  thermal  equilibrium  while 
still  in  motion.  The  temperature  at  which  the  base 
appears  to  have  stabilized  ranges  from  -C.6°C, 
when  the  ambient  temperature  was  -4.8°C,  to 
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G.2X,  when  the  ambient  temperature  was  --0.25CC. 
The  fact  that  temperatures  above  the  melting  tem¬ 
perature  of  the  snow  were  recorded  can  be  at¬ 
tributed  10  the  warm  air,  high  solar  input  and  a 
possible  offset  temperature  in  the  instrumentation. 
It  is  not  known  whether  thermocouple  3  in  Figure 
Cl  had  an  offset  temperature  because  it  was  de¬ 
stroyed  before  it  could  be  calibrated.  It  had  to  be 
destroyed  because  it  gave  erroneous  values  during 
later  tests.  These  plots  indicate  that  although  the 
speed  of  the  ski  was  still  increasing,  the  base 
reached  a  quasi-steady  condition  because  of  the 
generation  of  sufficient  frictional  heat  to  promote 
melting  over  the  portion  of  the  interface  that  was 
in  contact. 

This  observation  was  apparent  only  with  the 
Blitzski,  because  we  knew  the  location  of  the  ther¬ 
mocouples  in  the  Butz  ski  to  be  on  the  base.  In  the 
DH  ski,  the  thermocouples  were  located  nbov  e  t)  ie 
frictional  interface,  a  small  but  unknown  distance. 
This  small  distance  was  enough  to  severely  dampen 
their  thermal  response,  which  may  also  have  been 
affected  by  the  greater  thermal  cond  urti  vity  of  the 
DH  ski.  Nevertheless,  if  speeds  greater  than  18  m/ 
s  and  longer  runs  had  been  achieved  with  the  DH 
ski  at  similar  ambient  temperatures  (i.e.,  -5  to  CPC), 
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Figure  3-  Steady-state  vs  anwient  temperature  for  the 
Blitz  ski. 


similar  plots  showing  steady  temperatures  might 
have  been  observed.  Figure  3  shows  a  linear  rela¬ 
tionship  between  the  ambient  and  steady-state 
temperatures  as  calculated  from  the  data  given  in 
Figures  Cl  through  C4. 

The  discovery  that  the  base  of  a  slider  moving 
over  snow  can  attain  a  stead  y  temperature  supports 
the  idea  that  meltwater  at  the  frictional  interface  is 
responsible  for  the  low  friction  of  snow.  Clearly, 
McConica's  ideas  about  the  thermal  response  of 
skis  are  incorrect.  In  addition,  the  rfacial 
temperature  appears  to  be  independent  of  the 
speed  once  the  steady-state  temperature  has  been 
reached.  This  idea  is  illustrated  qualitatively  in 
Figure  4. 

Bowden  and  Hughes  (1939)  theorized,  and  it  is 
often  observed  that,  the  coefficient  or  friction  in¬ 
creases  with  decreasing  ambient  temperature.  The 
reason  for  this  is  directly  shown  by  the  lower 


t  teady-state  temperatures  that  we  found  at  lower 
ambient  temperatures.  Coulomb  or  dry  friction, 
which  has  a  higher  coefficient  than  the  lubricated 
friction,  is  relatively  more  important  at  lower  am¬ 
bient  temperatures.  Since  the  rate  of  heat  genera¬ 
tion  is  proportional  to  the  coefficient  of  friction,  it 
appears  that  the  greater  rise  in  interfacial  tem¬ 
perature  at  lower  ambient  temperatures  is  ascrib- 
able  to  more  heat  production  resulting  from  a 
higher  coefficient  of  friction.  This  effect  can  be  seen 
in  other  runs  that  did  not  reach  a  steady-state  tem¬ 
perature.  Comparing  Figures  Cl  7  and  Cl  8  shows 
thatforan  ambient  temperature  of -9.8°C,  the  tem¬ 
perature  rise  at  the  base  of  the  ski  was  5.3°C, 
whereas  for  an  ambient  temperature  of-6.7°C,  the 
temperature  rise  at  the  base  of  the  ski  was  only 
3.3°C.  Apart  from  the  change  in  ambient  tem¬ 
perature,  these  two  runs  were  nearly  identical. 

To  quantify  the  heat  flow  as  a  function  of  the 
ambient  temperature,  data  obtained  from  the  base 
of  the  skis  were  fitted  against  the  solution  for  the 
sudden  onset  of  a  constant  heat  flow.  From  Carslaw 
and  laeger  (1959,  p.  75) 

T(t)  =  T0A{%L\'/7  (8) 

ki\n  } 

where  TQ  is  the  initial  or  ambient  temperature. 
When  the  initial  temperature  rise  w  as  fit  ted  against 
t]/2  the  resulting  approximation  always  had  a 
close  fit.  The  values  of  the  heat  flow  constant,  or 
2iJs(&s/Jt)1/2/  k&,  are  shown  versus  ambient  tem¬ 
perature  in  Figure  5.  This  shows  that  the  heat  P  ow 
at  the  base  of  the  ski  increases  linearly  as  the  am¬ 
bient  temperature  decreases.  At  least  over  this  nar¬ 
row  range  of  temperatures,  the  friction  and  heat 


temperature  for  the  Blitz  ski. 


Ambient  Temperature  ( "C) 


Figure  5.  FIcat  flow  constant  vs  ambient  temperature 
for  the  Blitz  ski. 
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generation  increase  with  decreasing  ambient  tem¬ 
perature. 

Effects  of  load 

Fluctuation  of  load  at  constant  speed.  Figure  CIO 
illustraces  the  influence  of  load  cycling  on  heat 
generation.  The  thermal  response  resulting  from  a 
change  in  load  from  0.5  N  to  N  at  a  modest  speed 
was  surprisingly  large.  It  is  interesting  to  note  that 
starting  after  about  50  seconds,  a  curve  drawn 
through  the  peaks  would  approximate  f1/2.  That 
is,  although  the  load  oscillates  between  the  N/2 
and  N,  the  resulting  temperature  response  does 
not  oscillate  back  and  forth  between  the  maximum 
temperature  and  the  snow  temperature  because 
heat  accumulates  in  the  ski  base. 

This  effect  is  also  shown  in  FigureC19,  in  which 
the  ski  accelerates  downnill,  stops  and  accelerates 
downhill  again  in  three  cycles.  The  time  between 
stopping  and  starting  again  was  insufficient  for 
the  base  to  reach  a  state  of  therm;  1  equilibrium. 
Therefore,  when  heat  accumulates  as  the  ski  accel¬ 
erates  again,  the  ski  attains  a  higher  temperature 
than  it  did  during  the  previous  heating  episode. 

Different  constant  loads  during  descent.  Using  the 
iesponses  from  thermocouple  3  shown  in  Figures 
C8  and  C23,  as  well  as  some  of  the  data  not  included 
in  Appendix  C,  we  found  an  almost  linear  rela¬ 
tionship  between  the  load  applied  and  the  temper¬ 
ature  rise.  In  each  case  the  ambient  temperature 
was  approximately  -9.5°C.  This  is  in  accordance 
with  the  eq  1,  which  suggests  that  the  rate  of  heat 
generation  is  proportional  to  the  applied  load. 
This  relationship  is  illustrated  in  Figure6,  in  which 
the  no-load  data  point  is  assumed  to  produce  no 
temperature  rise. 

Turning.  Figures  C5  and  C6  show  the  tempera¬ 
ture  rises  associr.ted  with  the  turning  of  the  ski. 


Figure  C5  shows  the  response  of  a  basal  thermo¬ 
couple  during  a  downhill  run  on  a  long,  gentle 
slope.  As  is  quite  evident,  the  steady-state  temper¬ 
ature  is  approached,  but  is  never  reached  because 
of  the  weighting  and  unweighting  during  turning. 
Figure  C6  shows  the  therma’  response  of  two  basal 
thermocouples  that  are  located  on  either  side  of 
the  base  of  the  Blitzski.  The  fast,  smooth  tumspro- 
duced  a  significant  thermal  response  from  both 
sides  of  theski.  It  is  interesting  to  note  two  features 
in  this  figure.  First,  the  alternate  weighting  and 
unweighting  on  each  side  of  the  ski  during  turns 
produced  a  180°  phase  lag  between  the  two  basal 
thermocouples.  Second,  the  average  thermal  re¬ 
sponse  of  therm  ocouple  1  is  greater  than  that  of  5. 
This  can  be  attributed  to  the  inside  location  of  ther¬ 
mocouple  1  and  the  outside  location  of  thermo 
couple5.  In  downhillskiing,  more  weight  is  applied 
to  the  downhill  ski  during  a  turn.  In  the  case  of  the 
Blitz  ski,  which  was  placed  on  the  right  foot  for  all 
tests,  the  left  inside  edge  would  experience  more 
loading  when  negotiating  a  left-hand  turn  than  the 
outside  edge  would  experience  during  the  same 
turn.  Therefore,  thermocouple  1  was  more  heavily 
weighted  than  thermocouple  5  and  shows  a  greater 
thermal  response. 

Loading  with  skis  parallel.  It  wassurprising  to  dis¬ 
cover  the  large  temperature  rise  across  the  ski 
when  the  load  appeared  to  be  evenly  distributed 
during  a  parallel  schuss.  Figures  Cll  and  C12 
show  this  effect,  which  is  caused  by  the  amount  of 
edging  that  a  skier  applies  when  in  the  parallel 
tuck  position.  This  nonuniform  distribution  of 
weight  can  be  attributed  to  the  technique  required 
to  maintain  stability  when  skiing  in  that  position. 
The  amount  of  edging  required  appears  to  be  a 
function  of  the  hardness  of  the  snow  as  illustrated 
in  Figure  7.  This  graph  indicates  that  the  tempera- 
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Figure  7.  Temperature  rise  vs  transverse  position  for  tivo 
snow  types. 
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ture  is  more  uniformly  distributed  across  the  ski 
when  the  snow  has  a  low  hardness.  This  is  appar¬ 
ently  attributable  to  the  more  even  distribution  of 
load  when  the  soft  snow  conforms  to  the  shape  of 
the  ski.  A  similar  effect  of  soft  snow  will  be  shown 
later  when  we  discuss  the  longitudinal  profiles. 

Effect s  of  speed.  This  parameter  was  by  far  the 
hardest  to  control  or  record,  largely  because  of  the 
variations  in  terrain,  lack  of  an  appropriate  meas¬ 
uring  device  and  supporting  personnel.  However, 
the  thermal  response  as  a  function  of  speed  is  de¬ 
scribed  in  a  qualitative  sense  in  Figure  C22.  The 
la»  ge  temperature  rise  of  approximately  4°C  at  the 
beginning  of  this  run  is  from  a  maximum  speed  of 
approximately  6  m/s.  At  60  seconds,  a  smaller 
temperature  rise  of  approximately  2°C  is  due  to 
the  reduced  speed  on  the  Poma  lift. 

Effects  of  snow  type.  To  further  determine  the 
effects  of  snow  hardness  on  the  heat  production 
rate,  the  longitudinal  profile  was  analyzed.  This 
yielded  some  interesting  results,  as  shown  in  Fig¬ 
ures  C13  through  C16.  Figures  03  and  Cl4depic. 
the  longitudinal  thermal  response  at  the  base  of 
the  ski  on  hard-packed  snow.  Close  inspection  of 
the  thermocouple  response  associated  with  each 
curve  shows  that  in  icy  conditions,  most  of  the  heat 
generation  occurs  in  the  vicinity  of  the  center  of  the 
ski  (i.e.,  directly  under  the  skier).  In  the  case  of  soft 
snow  (Fig.  C15  and  C16),  the  weight  is  distributed 
far  more  evenly,  probably  because  the  soft  snow 
conformed  to  the  shape  of  theski.  The  temperature 
rise  as  a  function  of  longitudinal  position  for  hard 
and  soft  snow  is  illustrated  in  Figure  8. 

The  appn  ent  heat  flow  distribution  in  the  lon¬ 
gitudinal  direction  caused  us  to  nose  some  ques¬ 
tions.  Is  the  temperature  increase  at  the  rear  of  the 
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Figure  8 ■  Temperature  rise  vs  longitudinal  position  for 
different  snow  types. 


ski  or.  soft  snow 3  ttributable  to  an  accumulation  of 
heat  along  the  ski  or  is  it  ascribable  to  the  load  dis¬ 
tribution?  An  experiment  was  conducted  in  which 
the  static  load  distribution  for  the  DH  ski  in  the 
longitudinal  direction  was  measured.  To  determine 
the  pressure  distribution  along  the  length  of  the 
ski,  six  plastic  disks  were  machined  with  cross- 
sectional  areas  oil  6.4  mm2.  This  value  was  chosen 
so  that  the  pressure  applied  on  the  disks  would  be 
large  enough  to  be  measurable  on  the  Fuji  pressure 
film  that  was  available  to  us.  The  disks  were  equally 
spaced  along  thebase  of  theDHskiand  the  pressure 
film  was  placed  between  the  disks  and  a  hard,  Pat 
surface.  The  weight  of  the  skier  was  then  applied 
for  a  couple  of  minutes  until  the  film  responded. 
This  procedure  was  repeated  to  give  the  average 
pressure  of  the  two  tests  at  each  node  as  illustrated 
in  Figure  9. 

This  graph  indicates  that  there  is  more  pressure 
at  the  rear  than  the  front  of  the  ski,  so  more  heat 
would  be  generated  at  this  location.  In  addition, 
this  graph  indicates  that  the  largest  load  is  applied 
directly  under  the  skier,  which  explains  why,  on 
hard  snow  at  least,  a  large  thermal  response  was 
measured  under  the  skier.  It  appears  that  the  larg¬ 
est  thermal  response  occurs  at  Hip  rear  of  thp  ski  in 
soft  snow  conditions,  and  this  nuph  t  be  attributed 
to  the  accumulation  of  heat  along  the  length  of  the 
ski  as  well  as  the  increased  pressure  at  this  location 
under  soft  snow  conditions. 

Effects  of  thermal  conductivity.  Figure  C/  shows 
die  response  at  two  positions  in  the  Blitz  ski.  Ther¬ 
mocouple  1  was  on  the  base  and  responded  as  we 
have  seen  before.  Thermocouple  2  was  3.8  mm 
above  the  base  and  responded  much  less  and 
much  later.  The  difference  between  the  two  can  be 
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Figure  10.  Temperature  rise  vs  vertical  position  in  the 
Blitz  ski  at  —9.5°C. 


attributed  in  part  to  the  poorly  conducting  plastic 
of  che  Blitz  ski.  In  Figure  C8  the  responses  of  three 
thermocouples  cn  the  base  and  two  above  the  base 
are  shown.  Thermocouple  4,  which  was  9.7  mm 
above  the  base,  hardly  responded  to  the  hear  gen¬ 
eration.  It  is  interesting  to  note  that  the  temperature 
at  this  location  falls  off  at  the  initiation  of  motion. 
This  effect  can  be  attributed  to  convective  air  cool¬ 
ing,  implying  that  the  top  of  the  ski  was  not  at  the 
ambient  temperature,  even  though  it  appeared  to 
be  in  thermal  equilibrium.  This  convective  effect 
could  have  occurred  because  the  top  of  the  ski  re  ¬ 
ceived  heat  either  from  the  sun  or  from  the  skier's 
boot.  Figure  10  shows  the  dampening  of  the  thermal 
response  with  height  of  thermocouples  2, 4  and  5 
of  Figure  C8. 

Figure  C9  shows  the  vertical  response  of  the  DH 
ski  with  thermocouples  located  at  approximately 
the  same  vertical  positions  as  in  the  Blitz  ski.  Note 
that  tlu  relative  thermal  response  of  the 
thermocouples  away  from  the  base  (2a22  and  2a32) 
are  much  greater  than  the  response  in  the  Blitz  ski. 
This  is  because  thermocouple  2a22  is  situated  on 
top  of  the  aluminum  plate,  as  shown  in  Appendix 
B.  The  high  thermal  conductivity  of  this  plate 
increases  the  thermal  response  so  much  that  a 
thermocouple  off  of  the  base  in  the  DH  ski  (2a22) 
has  about  the  same  fractional  response  as  a 
thermocouple  on  the  base  (5)  of  the  Blitz  ski,  as 
shown  in  Figure  11. 

A  comparison  of  Figures  C20  and  Cll  shows 
the  variation  of  response  time  with  thermal  con¬ 
ductivity  during  the  same  ran.  (The  ambient  tem¬ 
perature  differs  in  these  figures  because  of  a  cali¬ 
bration  error.)  In  Figure  C20,  the  two  thermo¬ 
couples  on  the  steel  edge  (3c  10  and  4cl0) '  m  to 
respond  almost  immediately,  whereas  tin.  base 
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Figure  1 1 .  Temperature  rise  vs  vertical  position  in  the 
DH  ski  at  -2 .8  °C. 

thermocouple  (5c02)  responded  after4  or  5  seconds. 
In  Figure  Cll,  which  shows  only  base  thermo¬ 
couples,  the  response  begins  more  quickly  but 
only  2b01  responded  as  much  as  3cl0.  Response  is 
delayed  not  only  because  of  the  effect  of  thermal 
conductivity  but  alco  because  the  base  thermo¬ 
couples  in  the  DH  ski  are  a  small  but  finite  distance 
away  from  the  base. 

Additional  observations 

Variation  of  temperature  with  altitude.  FigureC17 
shows  the  change  in  temperature  that  occurs  with 
altitude.  In  this  case,  the  snow  temperature  is  ap¬ 
proximately  1 .3°C  lower  at  the  top  of  the  run  than 
it  is  at  the  bottom.  This  observation  indicates  that 
a  small  component  of  the  temperarure  rise  that 
takes  place  during  a  descent  may  be  caused  by  this 
altitude  effect. 

Real  location  of  basal  thermocouples  in  Blitz  ski.  Fig¬ 
ure  C8  shows  a  very  large  difference  in  thermal  re¬ 
sponse  among  the  three  basal  thermoc  -pies  (1  on 
the  lef‘  side,  2  in  the  middle  and  5  on  the  right  side 
of  the  Blitz  ski.)  This  might  be  attributable  to  the 
unevenload  distribution  in  the  transversedirection 
during  parallel  skiing  as  described  earlier.  This 
cannot  be  the  case,  however,  because  this  explan¬ 
ation  would  require  that  thermocouple  1  have  the 
strongest  response  and  thermocouple  5  have  the 
weakest.  The  reason  for  the  different  thermal  re¬ 
sponses  is  that  these  thermocouples  are  not  locat¬ 
ed  at  precisely  the  same  heights  above  the  base  of 
the  ski,  an  effect  that  is  especially  important  in  the 
nou-conductive  plastic.  Figure  C21  shows  thesame 
response  but  with  the  time  axis  expanded.  In¬ 
spection  of  the  response  time  suggests  that  ther- 
n .orouple  3  is  closest  to  the  base,  followed  by  5  and 
1.  lo  correct  this  problem,  new  thermocouples 


were  placed  in  the  base  of  the  Blitz  ski  and  a  typical 
response  of  these  is  shown  in  Figure  C22. 

Smoothing  of  data 

Figure  C24  shows  the  response  of  a  basal  thermo¬ 
couple  as  actually  recorded.  The  fluctuations  in 
thermal  response  during  the  run  may  be  caused  by 
electronic  noise  as  wellas  physical  effects,  including 
frequent  variations  in  load  due  to  the  topography 
of  the  ski  slope.  These  undesirable  effects  were  re¬ 
duced  using  a  smoothing  routine  with  a  cut-off 
frequency  of  0.5  Hz,  which  allowed  earlier  com¬ 
parisons  b  .'tween  files.  A  smoothed  version  of  Fig¬ 
ure  C24  is  shown  in  Figure  C23  for  compa.  son. 

FINITE  ELEMENT  MODEL 

Computational  procedure 

Purpose  of  computer  simulation  of 
frictional  heat  floiu  in  skis 

As  a  second  phase  to  this  project,  we  used  a 
computer  model  to  simulate  the  flow  of  frictional 
heat  into  a  downhill  racing  ski  to  try  to  quantify  the 
significance  of  thermal  conductivity  and  material 
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The  computer  model  made  use  of  the  finite  element, 
method  to  determine  the  temperaturedistribution 
in  a  rectangle  that  represented  the  cross  section  of 
tne  ski.  The  temperature  fields  generated  by  this 
model  were  compared  with  the  temperatures  re¬ 
corded  in  the  field  to  see  if  the  model  could  be  a 
use  :1  tool  to  predict  the  heat  flow  characteristics 
of  various  ski  structures  and  ambient  conditions. 
This  would  eliminate  the  expensive  and  tedious 
task  of  building  and  testing  sliders  of  all  types  to 
determine  their  thermal  response. 

Finite  element  software  ard  parameters 

The  THERMAP  (Thermal  Analysis  Program) 
software  available  on  the  NORTHS'! AR  computer 
network  at  Dartmouth  College  (Glovsky  1982) 
was  chosen  for  this  computer  simulation  because 
of  its  flexibility  in  handling  a  wide  variety  of  ther¬ 
mal  phenomena,  includ  ing  conduction,  convection, 
internal  heat  generation,  surface  heat  flux  and 
phase  chare  \  in  addition,  the  program  is  capable 
of  calculate  ,a  the  transient  and  steady-state  solu¬ 
tions  of  the  temperature  distribution,  both  of  which 
are  of  great  interest  in  frictional  heating.  Many  of 
the  input  parameters  were  based  on  observations 
and  measurements  from  the  field  study. 

To  use  the  finite  element  method,  a  mesh  had  to 


be  created,  which  involved  breaking  up  the  cross 
section  of  the  DH  ski  into  a  number  of  small 
elements.  Making  use  of  symmetry,  we  modeled 
only  one  half  of  the  wid  th  of  the  ski  with  a  no-heat- 
flow  boundary  condition  imposed  on  the  right 
side  of  the  mesh  at  the  center-line  of  the  ski.  To 
determine  the  element  dimensions  and  shapes,  a 
section  was  cut  from  an  old  DH  ski  and  the  material 
dimensions  and  configuration  were  measured. 
Furthermore,  the  size  of  the  individual  elements 
was  based  on  thee:'.  >ected  heat  flow  in  the  medium. 
The  elements  were  small  in  locations  where  the 
temperature  gradients  were  large.  In  the  case  of  a 
ski  sliding  over  snow,  the  elements  were  smallest 
at  the  base  of  the  ski,  where  the  heat  is  generated. 
We  incorporated  185  triangular  elements  into  the 
mesh.  A  larger  number  of  elements  could  have 
been  used,  but  for  this  application,  the  additional 
resolution  in  the  mesh  would  have  required 
extremely  long  computer  runs  The  185  elements 
were  created  from  the  1 1 4  nodes  that  were  specified; 
the  node  coordinates  are  the  locations  at  which  the 
resulting  temperatures  are  calculated.  A  schematic 
of  the  mesh  chosen  for  the  frictional  heat  flow  in 
the  DH  ski  is  illustrated  in  Figure  12. 

For  the  first  17  seconds  of  each  run,  a  heat  flux 
boundary  condition  resulting  from  frictional 
heating  at  the  base  of  the  ski  was  specified.  The 
heat  generation  was  determined  using  eq  1  with 
half  of  the  skier's  weight  assigned  to  each  ski.  The 
coefficient  of  friction  was  taken  as  0.06,  the  speed 
was  estimated  as  18  m/s,  and  the  area  was 
calculated  as  0.16  m2  by  measuring  the  length  and 
width  of  the  apparent  contact  area  of  the  DH  ski. 
This  speed  is  close  to  the  upper  limit  of  that 
encountered  with  the  DH  ski  but  is  less  than  racing 
speeds.  Actual  contact  only  occurs  over  a  small 
fraction  (-4%)  of  thebase  of  theski.  Using  this  area 
would  increase  the  heat  flux  by  a  factor  of  25  but 
would  be  applied  to  only  4%  of  the  base  of  the 
mesh  in  the  model .  Assuming  that  the  snow  crys  tals 
come  in  contact  with  all  of  the  ski  at  some  time  and 
that  these  contacts  vary  randomly  and  frequently, 
we  deemed  it  reasonable  to  apply  the  heat  flux 
across  the  entire  base  of  the  ski  with  a  value 
corresponding  to  100%  contact  area. 

All  of  this  heat  flows  into  either  the  ski  or  the 
snow  since  molting  has  not  yet  begun.  Using  this 
balance,  eq  8  for  both  the  snow  and  the  ski,  and 
noting  that  the  surface  of  the  snow  and  the  ski  are 
at  the  same  temperature,  we  can  easily  show  that 
the  heat  flow  ratio  (R)  is  given  by 

r  =r(a),/2-  (9) 
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Figure  1 2.  Finite  element  mesh  with  mule  num¬ 
bers.  Dimensions  are  in  millimeters.  The  right- 
hand  edge  is  the  cciiterliiu  ( the  ski  and  the 
left-hand  edge  is  the  side  of  the  ski. 


Using  the  values  fc.  =  C.377  vV/ m  K,  L  =  0.552  VV / 
m  K,  Ks  =  10~7  m2/ s,  and  k.  =  1.3  rn2/s  yielded  a 
value  of  0.438  foi  R.  Therefore,  the  rate  of  heat  flux 
into  tne  ski  was  taken  as 

Ts  -  ~  WN  .  (10) 

This  was  the  value  used  in  all  cases  simulating  fric¬ 
tional  heating. 

As  an  alternative,  a  fixed  temperature  condition 
could  be  imposed  on  any  number  of  nodes  within 
the  mesh.  From  observations  made  in  the  field 
study,  wedetermined  that  once  theski  basereached 
a  steady  temperature,  this  steady-state  condition 
was  maintained  for  the  duration  of  the  run;  its 
value  in  most  runs  was  above  -1°C.  To  simulate 
this  temperature  condition  in  the  model,  the  pro¬ 
gram  was  stopped  once  the  base  had  attained  tem¬ 
peratures  within  this  range.  A  fixed  temperature 
was  then  imposed  on  the  basal  nodal  points,  and 
the  simulation  continued. 

In  the  case  where  a  rise  and  subsequent  fall  i” 
the  vertical  temperature  profile  was  determmed,  it 
was  possible  to  simulate  the  deceleration  of  the  ski 
using  appropriate  values  for  basal  temperature. 


This  was  achieved  by  decrement  n.g  the  tempera¬ 
ture  of  the  basal  nodes  in  a  way  that  approximately 
replicated  the  temperature  decreases  observed  in 
the  field  data  during  deceleration  of  the  ski-  Once 
again,  we  did  this  by  stopping  the  program  to  ad¬ 
just  the  fixed  temperatures  of  the  basal  nodes  to 
simulate  the  cooling  effect  observed  in  the  field. 

Internal  heat  generation  from  flexing  could  have 
been  included  but,  from  the  data  gathered  in  the 
field,  there  is  no  indication  that  heat  was  being 
generated  at  any  location  other  than  the  base  of  the 
ski.  Tire  effects  of  convection  ascribable  to  motion 
occur  along  those  boundaries  that  were  exposed  to 
the  air,  the  left  and  top  side  of  the  mesh.  The  con¬ 
vection  coefficient  was  determined  using  the  equa¬ 
tion  (Kays  1960) 

H  =  7.25  u4/5  (ID 

where  the  speed  was  taken  as  18  m/s  (H  being 
given  in  W/m2  °C).  This  value  was  used  for  all 
simulations.  Values  of  the  properties  of  the 
materials  used  in  the  mesh,  including  the  thermal 
conductivity  and  diffusivity,  are  gUen  in  Table  1. 

Case  studies  of  heat  flow  into  the  DH  ski 

We  used  four  case  studies  to  investigate  the 
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Table  1.  Thermal  properties  of  materials 


in  DH  ski. 
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effect  of  thermal  conductivity  on  the  retention  of 
heat  at  the  base  of  a  slider.  For  each  case,  both  the 
transient  and  steady-staU  solutions  for  the 
temperaturodistributions  were  determined.  The  .e 
cases  are  as  follows. 

Gtse  ] .  The  temperature  distribution  was  calcu¬ 
lated  for  the  actual  RossignolDH  ski,  incorporating 
both  the  aluminum  plate  near  the  base  and  the 
steel  edge.  This  case  gave  the  heat  flow  character¬ 
istics  of  this  popular  downhill  racing  ski.  This  con¬ 
figuration  is  shown  in  Appendix  13. 

Case  2.  The  temperature  distribution  was  calcu¬ 
lated  for  the  OH  ski  with  the  steel  edge  replaced  by 
a  ceramic  edge.  The  removal  of  the  highly  conduc¬ 
tive  steel  edges  might  allow  greater  melting  and 
give  lower  friction  at  low  temperatures.  This  could 
be  achieved  by  substituting  a  high  strength,  poorly 
conducting  material,  such  as  sintered  alumina,  for 
the  steel. 

Cased.  The  temperature  distribution  was  calcu¬ 
lated  with  the  steel  edge  replaced  with  the  ceramic 
edge  and  the  lower  aluminum  plate  replaced  with 
a  layerof  polymer.  The  vibration-dampening  alum¬ 
inum  layer  used  in  racing  skis  was  thought  to  be  a 
significant  heat  flow  conduit  because  of  its  very 
high  thermal  conductivity,  geometry  and  location 
in  the  ski.  We  hoped  that  the  replacement  of  this 
aluminum  layer  with  a  material  of  much  lower 
thermal  conductivity  but  comparable  stiffness, 
such  as  polyoxymethylene,  would  reduce  the  heat 
flow  along  this  layer,  resulting  in  greater  retention 
of  heat  at  the  base. 

Case  4.  The  temperature  distribution  was  calcu¬ 
lated  for  a  ski  with  the  aluminum  layer  replaced 
with  the  polyoxymethylene  layer  but  with  the 
steel  eage  reinserted.  This  material  configuration 
was  analyzed  to  determine  if  the  two  metal  com¬ 
ponents  interacted  to  enhance  the  heat  removal. 

A  list  of  the  thermal  properties  of  materials  in 
the  actual  and  altered  DH  ski  is  given  in  Table  1. 


LimitutionstUhl  problems  encountered  williTHERMAP 

In  the  transient  solutions,  there  was  a  small 
amount  of  heat  flow  across  the  centerline  of  the  ski, 
which  was  designated  as  a  no-heat-flow  boundary. 
This  could  have  been  caused  by  the  selection  of 
mesh  size  and  element  configuration.  However, 
the  effects  of  this  problem  were  small  and  we 
considered  them  insignificant  because  our  interest 
was  primarily  in  the  heat  flow  along  the  base  and 
out  the  sides  of  the  ski. 

In  trying  to  simulate  the  temperature  response 
at  tire  base  of  tire  sk,  as  observed  in  the  field  study, 
we  used  a  heat  flux  boundary  condition  to  start 
the  computation,  as  described  above.  Its  use  was 
limited,  however,  because  the  temperature  at  the 
base  was  limited  to  the  observed  steady-state 
temperature.  Thus,  the  boundai  y  condition  at  the 
base  was  adjusted  according  to  the  calculated 
temperatures.  This  required  frequent  pauses  in 
the  execution  of  THERM  AP  to  adjust  the  heat  flux 
or  impose  fixed  temperatures  on  specific  nodal 
points.  This  procedure  is  somewhat  arbitrary  but 
is  not  unrealistic  since  the  friction  and  the  heat 
production  vary  with  time  after  the  start  of  motion 
in  a  real  ski.  It  would  have  been  convenient  to  set 
these  parameters  in  an  additional  file  that  would 
have  been  canahleofr.  iinstinr*  theheat  flux  values 
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and  the  convection  coefficients  at  the  specitied 
times  and  locations,  but  the  actual  heat  flux  could 
not  be  determined  from  data  collected  in  the  field. 
Apart  from  this  limitation,  the  THERM  AP  software 
package  was  an  appropriate  tool  in  determining 
the  temperature  response  resulting  from  frictional 
heatu.g. 

Results  and  discussion 

The  finite  element  heat  flow  analysis  was  done 
to  examine  the  significance  of  thermal  characteris¬ 
tics  to  the  heat  flow  propagation  in  a  slider.  The  re¬ 
sults  presented  and  discussed  in  this  section  yield 
some  interesting  insights  into  the  heat  flow,  and 
thus  we  hope  that  this  model  can  be  used  to  predict 
the  thermal  and  frictional  behavior  of  sliders. 

As  noted  above,  for  about  the  first  17  seconds 
the  heat  flow  parameters  assigned  were  based  on 
continuous  heat  generation  calculated  from  eq  10 
and  the  area  of  the  ski.  This  initial  period  of  constant 
heat  generation  does  not  take  into  account  the  ef¬ 
fects  of  weight  or  speed  change  that  occurred  in 
the  field  study  because  of  the  difficulty  of  adjusting 
the  prescribed  temperatures  at  nearly  every  time 
increment.  When  the  base  reached  a  temper.'  are 
that  was  comparable  to  the  steady-state  tempera¬ 
ture  measured  in  the  field,  that  temperature  was 
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fixed  to  the  basal  nodes  and  this  condition  was 
maintained  until  the  entire  ski  was  in  a  state  of 
thermal  equilibrium.  Since  the  replication  of  a 
common  ski  run  was  desired,  the  basal  tempera¬ 
tures  were  also  modified  to  simulate  the  cooling 
experienced  at  th e  base  as  the  ski  deed  ated  and 
stopped.  This  simulation  allowed  a  qualitative 
comparison  between  the  vertical  thermocouple 
response  in  the  model  and  the  test  ski. 

Transient  thermal  response  of  the 
actual  DH  ski  (case  1  j 

Figures  Dl  through  DIO  illustrate  the  transient 
heat  propagation  through  the  ski  until  a  steady- 
state  condition  has  been  achieved.  These  plots  of 
isotherms  show  how  the  temperature  distribution 
varies  with  time,  illustrating  the  long  time  required 
before  the  entire  ski  is  in  a  complete  state  of  thermal 
equilibrium.  To  show  this,  the  temperature  of 
node  46  (see  Fig.  12)  in  easel  was  plotted  asa  func¬ 
tion  of  time  and  is  shown  in  Figure  1 3.  This  suggests 
that  ‘he  time  required  for  the  entire  ski  to  achieve 
complete  steady-state  far  exceeds  the  time  of  a 
standard  downhill  competition  run.  Except  at  the 
base  of  the  ski,  the  id  ja  of  steady  state  may  not  be 
applicable  to  skis  because  of  the  frequent  changes 
iu  speed,  direction  and  load  experienced  under 
nearly  all  skiing  conditions. 

V  tical  thermal  response 

We  did  this  computer  run  to  determine  the 
thermal  characteristics  of  nodes  in  the  ski  that  ap¬ 
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Figure  13.  Temperature  riseai  node  46  vs  time  for  case  1 . 


proximately  represent  the  measured  vertical  tem¬ 
perature  profile  in  the  DH  ski.  It  should  be  noted 
that  the  computer  simulation  represents  an  ideal 
downhill  run  consisting  of  no  intermittent  heating, 
uniform  heating  across  the  base  of  the  ski,  and  a 
run  time  representing  the  time  of  a  normal  downhill 
competition  run.  Because  of  these  assumed  condi¬ 
tions,  only  a  qualitative  comparison  can  be  made 
between  model  and  field  results.  The  response  of 
the  computer  simulation,  shown  in  Figure  14,  is 
typical  of  the  amplitude-dampened,  phase-lagged 
behavior  observed  in  the  real  ski  and  shown  in 
Figure  C9. 

Thermocouple  na02  reached  a  plateau  at  a  tem¬ 
perature  of  approximately  -1°C,  thus  simulating 
what  was  observed  in  the  field  data.  The  thermo¬ 
couples  furthest  away  from  the  base  of  the  ski 
(na32  and  na42)  cool  off  before  experiencing  the 
frictional  heat  propagating  up  from  the  base.  In  the 
real  ski,  this  was  attributed  to  the  effects  of  convec¬ 
tion.  However,  in  the  model  this  can  only  be  attrib¬ 
uted  to  computational  inaccuracy. 

Heat  flux  at  the  base 

The  heat  flux  perpendicular  to  the  base  was 
computed  at  four  places  along  thebase  (see  Fig.  15) 
because  that  is  where  the  critical  heat  removal 
takes  place.  It  is  clear  a  orn  tl  iese  plots  that  the  local 
heat  flux  is  strongly  affected  by  material  compo¬ 
sition  and  is  very  dependent  on  the  location  along 
the  base.  As  expected,  heat  flux  increases  with  in¬ 
creasing  thermal  conductivity,  but  the  effect  is  not 
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Figure  14.  Computed  vertical  temperature  profile  for 
nodes  110,  76,  58,  40  (equivalent  to  thermocouples 
na02,  na22,  na32,and  na42  respectively). 
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Figure  15.  Heat  flux  versus  time. 


uniform  throughout  the  complex  geometry  of  this 
ski.  Case  1  incorporated  an  aluminum  plate  with  a 
steel  edge;  case  2  incorporated  an  aluminum  plate 
with  a  ceramic  edge;  case  3  incorporated  a  polymer 
plate  with  a  ceramic  edge,  and  case4  incorporated 
a  polymer  plate  with  a  steel  edge.  The  locations  of 
the  nodes  in  this  heat  flux  analysis  are  shown  in 
Figure  12. 

Figure  15a  shows  the  heat  flux  as  a  function  of 
time  between  basal  nodes  100  and  91,  which  are 
located  approximately  3  mm  from  the  edge  of  the 
ski.  The  heat  flows  in  cases  1  and  4  with  the  steel 
edge  are  about  four  times  greater  than  in  cas' s  2 
and  3  with  the  ceramic  edge,  at  least  at  a  time  of  227 
seconds.  This  suggests  that  the  steel  edge  dominates 
the  heat  flow  pattern  this  close  to  the  edge.  The 
retention  of  heat  near  the  side  of  the  ski  base  is  in¬ 
creased  significantly  when  a  material  of  low 
thermal  conductivity  is  used  in  place  of  steel. 

The  effects  of  the  steel  edge  and  the  aluminum 
plate  on  the  heat  flux  at  thebase  of  theski  are  about 
equal  at  a  distance  of  10  mm  from  the  edge.  In  the 


case  of  basal  nodes  101  and  92  (Fig.  15b),  both  the 
steel  edge  and  the  aluminum  plate  are  needed  to 
make  a  significant  difference.  Clearly,  the  heat  flux 
becomes  more  dependent  on  the  aluminum  plate 
than  on  the  steel  edge  as  the  distance  away  from 
the  edge  increases.  The  heat  flow  between  nodes 
1 1 0  and  107,  located  25  mm  away  from  the  edge,  is 
shown  in  Figure  15c.  The  heat  fluxes  at  227 seconds 
for  cases  3  and  4  are  very  similar  to  but  slightly  less 
than  that  of  case  2  (aluminum  plate,  ceramic  edge). 
The  aluminum  plate  really  dominates  other  effects 
at  37  mm  away  from  the  edge  as  shown  in  Figure 
15d  for  nodes  112  and  114.  Cases  1  and  2  with  the 
aluminum  plate  provide  the  greatest  amount  of 
heat  flow  and  would  reduce  the  amount  of  heat  at 
the  base  of  the  ski  in  its  midsection  by  50%  or  more. 

Steady-state  temperature  distribution 

Figure  16  illustrates  the  steady-state  tempera¬ 
ture  distribution  of  cases  1  through  4.  Cases  1  and 
2  (Fig.  16a  and  b)  with  the  aluminum  plate  show 
that  a  large  amount  of  heat  is  conducted  along  the 
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aluminum  plate  to  the  outside  wall  at  the  ski, 
whereas  in  cases  3  and  4  (Fig.  16c  and  d)  the  hori¬ 
zontal  component  of  heat  flow  in  the  polymer 
plate  is  significantly  reduced.  1  his  suggests  that, 
when  the  entire  ski  is  in  an  overall  state  of  thermal 
equilibrium,  a  considerable  amount  of  heat  is  lost 
through  the  side  of  the  actual  ski.  Figures  16c  and 
d  show  that  the  most  of  the  heat  flows  vertically  ip 
into  the  ski  and  that  the  presence  of  the  polymer  re¬ 
duces  the  total  heat  flow,  although  there  is  still  a 
significant  corner  effect  from  the  steel  edge  in 
Figure  16d. 

In  an  effort  to  quantify  this,  we  measured  the 
perpendicular  distances  between  adjacent 
isotherms  for  each  case  to  determine  the  magnitude 
of  the  heat  flow  in  the  plate  and  the  surrounding 
material.  Results  indicate  that  for  case  1,  the  heat 
flux  along  thealumnuim  plate  is  approximately  40 
times  as  great  as  the  heat  flow  into  the  wood  just 
above  the  aluminum  plate.  In  case  2,  the  replace¬ 
ment  of  the  steel  with  ceramic  reduces  this  ratio  to 
approximately  1 5.  This  red uction  occurs  in  spite  of 
the  downward  heat  flow  in  the  vicinity  of  the  steel 


edge.  In  cases  3  and  4  this  ratio  is  reduced  to  ap¬ 
proximately  1,  thus  showing  a  greatly  reduced 
heat  loss  out  of  the  sides  of  the  ski.  This  analysis 
clearly  suggests  that  the  replacement  of  conductive 
materials  would  reduce  the  heat  flow  away  from 
the  base  and  could  greatly  affect  the  ski’s  friction 
at  lower  temperatures. 

The  heat  flow  in  the  wood  core  is  not  greatly  af¬ 
fected  by  the  change  in  material  composition,  per¬ 
haps  by  no  more  20%.  In  case  1,  so  much  heat  is 
conducted  along  the  aluminum  plate  that  heat 
flowsdownward  through  the  steel  edge.  However, 
this  may  be  a  feature  attributable  to  the  boundary 
conditions  that  were  imposed  on  the  model. 

A  qualitative  comparison  between  the 
field  study  and  the  finite  element  model 

A  qualitative  comparison  was  made  between 
the  results  of  the  field  study  and  the  results  of  the 
model  partly  to  test  the  validity  of  the  model. 

Tnw>icnt  rcsinvisc 

For  the  first  1 7  seconds  of  a  model  run,  a  constant 
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heat  flux  was  imposed  on  the  ski  base.  In  the  field 
study,  it  took  a  certain  amount  of  time  to  attain  this 
maximum  speed,  which  would  suggest  that  the 
amount  of  heat  generated  at  the  base  increased 
from  zero  to  its  maximum  value.  However,  if  we 
look  at  the  response  of  the  basal  thermocouple  in 
Figure  C2,  the  time  required  to  reach  steady  state 
is  smail  compared  to  the  time  for  which  the  base  is 
in  astate  of  thermal  equilibrium.  If  the  basal  steady- 
state  lasts  for  a  minute  or  more,  as  was  assumed  in 
the  model,  then  the  heat  flow  conditions  imposed 
in  the  model  closely  resemble  the  actual  field  pro¬ 
cess.  This  is  further  supported  by  the  times  it  takes 
to  attain  steady  state  in  the  model  and  in  FigureC2, 
which  are  within  a  factor  of  two  of  each  other. 

Vertical  temperature  profile 

For  this  comparison,  thecoordinates  of  the  nodes 
in  the  finite  element  mesh  were  selected  by  the  ap¬ 
parent  locations  of  the  thermocouples  in  the  DH 
ski.  Both  Figures  0)  and  l4show  the  phase-lagged, 
amplitude-dampened  response  of  thermocouples 
located  above  the  frictional  heat  source.  It  is  hard 
to  decide  just  how  closely  to  compare  thermo¬ 
couples  2a02  and  na02  because  the  times  of  fric¬ 
tional  heating  in  the  field  and  computer  runs  were 
quite  different.  VViiiie  the  run  in  Figure  C9  lusieu 
about  30  seconds,  the  simulated  run  depicted  in  Fi¬ 
gure  14  lasted  for  over  3  minutes.  Nevertheless, 
the  general  characteristics  of  these  two  profiles  are 
similar. 

Simplification s  of  the  computer  model 

From  the  data  collected  in  the  field,  we  assumed 
that  the  weight  of  a  downhill  skier  is  not  evenly 
distributed  across  the  width  of  the  ski,  but  is  ap¬ 
plied  more  to  the  inside  edge.  This  results  in  more 
heat  generation  nenrtheedgeof  theski,  which  was 
not  taken  into  account  in  the  computer  model.  The 
prevailing  philosophy  is  that  it  is  desirable  to  have 
the  ski  flat  on  the  snow  so  that  meltwater  can  be 
evenly  generated  across  the  entire  width  of  theski. 
Therefore,  heat  was  generated  evenly  across  the 
base  of  the  ski  in  the  finite  element  simulation 
However,  this  hindered  comparisons  of  basal  heat 
flow  between  the  model  and  the  field  data. 

Also,  heat  generation  in  the  model  was  not  in¬ 
termittent,  as  it  was  in  the  field  because  of  changes 
in  speed  and  load.  Because  of  the  randomness  of 
those  variations  in  real  ski  runs,  we  decided  that 
incorporating  these  effects  would  complicate  the 
analysis  too  much. 


CONCLUSIONS 

It  was  encouraging  to  have  obtained  such  a 
large  thermal  response  from  frictional  heating  at 
the  base  of  both  test  skis.  We  initially  thought  that 
the  response  of  frictional  heat  generation  might  be 
lost  in  the  noise  of  the  data  acquisition  system.  In¬ 
stead,  the  large  thermal  response  has  allowed  us  to 
characterize  the  nature  of  heat  generation  and  heat 
flow,  and  helps  us  understand  the  coefficient  of 
kinetic  friction. Clearly,  someof  McConica's  (1 950) 
ideas  of  snow  friction  have  been  disproved  and 
some  ideas  of  Bowden  and  Hughes  (1939)  have 
been  affirmed. 

At  the  outset  of  this  investigation,  we  hoped 
that  simple,  but  careful  tests  of  sliders  passing 
over  snow  could  provide  meaningful,  quantitative 
results  that  could  be  used  to  test  the  accuracy  of  the 
tempoiature  ri'-e*  predicted  by  theory.  The  theory 
might  then  be  more  useful  in  quantifying  friction 
under  known  ambient  conditions.  However,  be- 
causeof  the  highly  varying  conditions  that  existed 
in  the  ski  tests,  it  was  difficult  to  control  the  para¬ 
meters  a  fleeting  the  rateof  frictional  heating.  These 
included  the  speed  and  the  weight  acting  on  the 
slider.  Although  this  detracted  from  the  quantita¬ 
tive  i  csLihs,  tile  qurutativ o  analysts  gives  us  some 
insight  into  how  the  frictional  h  it  is  distributed 
and  supports  existing  ideas  abt  t  the  meltwater 
iayer. 

The  field  results  support  the  following  conclu¬ 
sions: 

1 .  The  formation  of  a  discontinuous  meltwater 
layer  is  suggested  by  the  steady-state  temperature 
reached  at  the  base  of  the  ski.  This  temperature 
was  consistently  close  to,  but  less  than,  the  melting 
temperature  of  the  snow.  There  appears  to  be  an 
upper  limit  to  the  basal  temperature,  which  is 
determined  by  some  combination  of  the  fractional 
area  in  contact  and  the  melting  temperature. 

2.  The  thermal  response  is  dependent  on  the 
ambient  temperature  of  the  snow.  This  suggests 
that  more  meltwater  is  present  at  higher  ambient 
temperatures,  thus  the  lower  friction. 

3.  The  heat  generated  by  friction  at  the  interface 
is  sensitive  to  the  ambient  snow  tempet  .ure,  the 
load  and  the  speed  of  the  slider. 

4.  The  distribution  of  frictional  heating  is  sensi¬ 
tive  to  the  style  of  skiing  and  the  hardness  of  the 
snow  surface.  The  heat  generation  is  more  uni¬ 
formly  distributed  if  the  snow  is  soft  and  can  con¬ 
form  to  the  shape  of  the  ski.  The  heating  is  con¬ 
centrated  under  the  skier  when  the  snow  is  hard. 
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5.  The  steady-state  temperature  ol  the  base  of  a 
ski  appears  to  be  independent  of  the  speed,  prob¬ 
ably  because  the  average  temperature  at  tire  base 
is  determined  by  the  melting  temperature  and  the 
contact  area,  even  if  the  heat  production  increases. 

6.  When  in  a  tuck  position  a  skier's  weight  is 
concentrated  on  the  inside  edge. This  result  shows 
that  the  ski  is  not  Hat  on  the  snow  even  when  the 
skier  perceiv  es  it  as  such.  Temperature  measure¬ 
ments  should  be  used  to  learn  how  to  correct  this 
uneven  distribution  so.  as  to  reduce  friction. 

7.  The  vertical  temperature  profile  exhibits  the 
charactensticpb  -e-lagged,  amplitude-dampened 
response  expected  when  heat  is  generated  at  one 
surface  only. 

8.  Temperature  measurements  are  a  powerful 
tool  for  examining  such  things  as  the  effectiveness 
of  waxes  and  slider  designs  and  should  become  a 
standard  tool  for  testing  the  behavior  of  skis  and 
other  sliders. 

The  finite  element  model  results  support  the 
following  conclusions: 

1 .  A  large  amount  of  heat  flow  occurs  out  the 
side  of  tire  DH  ski  because  it  has  the  higlrlv  con¬ 
ductive  aluminum  plate.  The  heat  flux  along  the 
aluminum  plate  is  about  40  times  greater  than  up 
through  the  wood.  Considering  that  the  aluminum 
plate  is  about  one-fortieth  as  thick  as  tire  wood  is 
wide,  about  the  same  amount  of  heat  flows  out  the 
aluminum  plate  as  flows  up  through  the  wood. 
Thus,  the  aluminum  plate  about  doubles  the  heat 
loss  from  the  ski.  The  replacement  of  this  plate 
with  a  material  of  lower  thermal  conductivity  but 
equal  stiffness,  such  as  polyoxyniethyleiro,  would 
inhibit  the  flow  of  heat  away  from  the  frictional 
interface  and  increase  the  supply  of  meltwater  at 
the  base  of  the  ski. 

2.  l  ire  heat  flow  out  the  sides  of  the  ski  is  in¬ 
creased  by  tire  steel  edge,  but  not  to  the  same  ex¬ 
tent  as  by  the  aluminum  plate.  However,  the  two 
complement  one  another  in  extracting  heat  from 
the  base  of  the  ski,  and  the  steel  edge  is  very  effec¬ 
tive  at  removing  heat  at  the  corner  itself. 

3.  lire  vertical  temperature  profile  exhibits  the 
characteristic  phase-lagged,  amplitude-dampened 
response  seen  in  the  tests  skis,  thus  showing,  that 
both  respond  to  heat  at  the  base. 

4.  The  finite  element  model  is  a  useful  tool  for 
predicting  the  thermal  response  of  sliders  and  can 
be  used  to  predict,  in  a  qualitative  way,  the  frictional 
characteristics  of  anything  that  slides  on  snow. 
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APPENDIX  A:  CR 10  DATA  ACQUISITION  PROGRAMS 


Key  strokes 

Comments 

Keystrokes  Comments 

Example  program  used  to  sample 

A 

9951 

thermocouples  in  the  Blitz  ski: 

A 

9999 

* 

1 

Access  program  table  #1 

A 

87 

A 

1 

Set  1  Hz  sampling  frequency 

A 

0 

A 

11 

Record  reference  junction 

A 

32 

temperature 

A 

86 

A 

1 

A 

76 

A 

1 

A 

14  Record  thermocouple  tempera¬ 

A 

3 

tures  connected  to  AM  32  board 

A 

1 

A 

1 

A 

1 

A 

1 

A 

0 

A 

2 

A 

14 

Record  thermocouple 

A 

1 

temperatures 

A 

1 

A 

5 

A 

2— 

A 

1 

A 

1 

A 

2 

A 

0 

A 

1 

A 

95 

A 

1 

A 

20 

A 

2 

A 

9990 

A 

A 

1 

n 

A 

A 

9999 

1.1  D _ J  l1 _ _ _ 1  ~  - 

A 

U 

86 

Move  all  values  to  final  memory 

r\ 

l'*  ANCCCUU  UlCiiilCICCJUpiC  ICilipCld” 

tures  connected  to  wiring  panei 

A 

10 

A 

4 

A 

70 

Place  all  values  in  storage  module 

A 

1 

memory 

A 

3 

A 

5 

A 

1 

A 

2 

A 

1 

A 

96 

A 

36 

A 

71 

A 

1 

* 

0 

Compile  and  run  program 

A 

0 

A 

86  Move  all  values  to  final  memory 

Example  program  to  sample 

A 

10 

thermocouples  in  the  DH  Ski: 

A 

70  Place  all  values  in  storage  module 

* 

A 

Access  input  number  address 

memory 

39 

A 

Set  number  of  inputs  to  39 

A 

39 

* 

1 

Access  program  table  #1 

A 

2 

A 

2 

Set  0.5  Hz  sampling  frequency 

A 

96 

A 

11 

Record  reference  junction 

A 

71 

temperature 

* 

0  Compile  and  run  program 

A 

1 

A 

1 

Keystrokes  used  to  manipulate  data: 

A 

3 

♦ 

6  Observe  real  time  measurements 

A 

1 

* 

7  Observe  measurements  sent  to 

A 

1 

storage  module  memory 

A 

0 

* 

9  Place  filemarksand/orerasedata  in 

A 

20 

storage  module  memory. 

APPENDIX  B:THERMOCOUPLE  AND  MATERIAL  CONFIGURATIONS  OF  THE  SKIS 


1.1  rn 


Figure  Bl.  Longitudinal 
position  of  all  thermocoup¬ 
les  in  the  Blitz  Ski. 


#1 ,  #2,  #3, 
#4,  #5 


0.5m 


O  HOMOGENEOUS  PLASTIC. 


Figure  B2.  Vertical  and  transverse  locations  of  ther¬ 
mocouples  for  test  runs  Fieldl  I  .dat  through  Field75.dat, 
FieldAl.dat  through  FieldA35.dat,  and  SKl-El.dat 
through  SKI-E9.dat  with  the  Blitz  Ski. 
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Figure  B3.  Vertical  and  trans¬ 
verse  locations  of  thermocouples 
for  test  runs  Ficld8l  .dat  through 
Field918.dat  with  the  Blitz  Ski. 
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Figure  B4.  Longitudinal  position  of  all  thermocouples  in  the  DH  Ski. 
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Figure  65.  Location  of  base  and  edge  thermocouples  ar 
the  front  and  rear  of  the  DH  ski. 
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Figure  B6. Transverse  profile  of  the  basal  thermocouples  located  in 
the  2bnn,  3bnn  and  4bnn  longitudinal  positions  in  the  DH  ski. 
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Figure  87.  Vertical  profile  of  the  thermocouple $  located  in  the  2ann,  3niui  mid 
4ann  longitudinal  positions  in  the  DH  ski. 
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Figure  BS.  Position  of  the  edge  thermocouples  located  in  the  2 enn,  3cnn 
and  4cnn  longitudinal  positions  in  the  DH  ski. 


APPENDIX  C:  TEMPERATURE  VERSUS  TIME  PLOTS  FROM  THE  FIELD  DATA 


Samples  of  the  temperature  measurements  are 
given  here.  All  of  the  curves  have  been  smoothed 
with  a  0.5-Hz  cut-off  frequency.  The  initial  tem¬ 
perature  at  the  base  of  the  ski  is  taken  as  the  ambi¬ 
ent  temperature.  The  test  run  and  snow  conditions 
are  described  briefly.  The  start  of  motion  is  indi¬ 


cated  by  A  and  the  stopping  point  by  B.  The  ther¬ 
mocouples  are  indicated  by  number;  the  location 
of  each  thermocouple  in  the  ski  is  shown  in  Appen¬ 
dix  B.  The  ski  can  be  identified  by  the  number  of 
the  thermocouple. 


Figure  Cl.  File  name:  Field  3A.dat.  Ami  lent 
temperature:  -0.25°C.  Maximum  speed:  =7  m/s. 
Run  description:  Parallel  down  short,  gentle 
incline  with  load  of  =N/2.  Snow  condition:  Soft, 
fresh  snow. 


'■me  [ s! 


Figure  C2.  File  name:  Field  5C.dat.  Ambient  ” 
temperature:  -4.8°C.  Maximum  speed:  =9  m/s.  ~ 
Run  description:  Parallel  down  short, gentle  incline  I 
with  load  of  =N/2  Snow  condition:  Hard,  old,  dry  f 
snow. 


Time  (S ) 


Figure C3.  File  name:  Ski-E5.dat.  Ambient  temperature: 
-3.4°C-  Maximum  speed:  =22  m/s.  Run  description: 
Parallel  down  short,  steep  incline  with  load  o/=3N/4. 
Snow  condition:  Soft,  fresh  snow. 
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Figure  C4.  File  name:  Ski-E9.dat.  Ambk.it  temper¬ 
ature:  -3.5°C.  Maximum  speed:  8  mis.  Run  de¬ 
scription:  Parallel  down  long,  gentle  incline  with 
load  of  =N/2.  Snow  condition:  Soft,  fresh  snow. 
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Figure  C5.  File  name:  Ski-E6.dat.  Ambient  tem¬ 
perature'  -4.75°C.  Maximum  speed:  -9  m/s.  Run 
description:  Simulation  of  a  downhill  run.  Time  of 
run  ~75  seconds.  Snow  condition:  Soft,  fresh 


~~  Figure  C6.  File  name:  Ski-E2.dat.  Ambient  tem¬ 
perature:  -3.2"C.  Maximum  speed:  =0  m/s.  Run 
description:  Wide, smooth  turns.  Snow  condition: 
Soft,  fresh  snow. 


Figure  C7.  File  name:  Ski-E3.dat.  Ambient  tem¬ 
perature:  -3.9°C.  Maximum  speed:  ~12  m/s.  Run 
description:  Parallel  down  short,  steep  incline  with 
load  u/=N/2.  Snow  condition:  Soft,  fresh  snow. 


Figure  CS.  File  name:  Field  410.dat.  Ambient  tempera¬ 
ture:  -9.5°C.  Maximum  speed:  =6  m/s.  Run  description: 
Parallel  down  abort,  gentle  incline  with  load  of  =N. 
Snow  condition:  Dry,  packed  natural  snow. 


Time  (s ) 


FigureC9.F He  name:  Field  C4D.dat.  Ambient  temperat¬ 
ure:  -1  .S°C.  Maximum  speed:  =9  m/s.  Run  description: 
Parallel  down  short,  gentle  incline  with  load  of  =N/2. 
Snow  condition:  Soft,  fresh  snow. 
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Figure  Cl 0.  File  name:  Field  S3.dat.  Ambient  tempera-  | 
ture:-l  ,S°C.  Maximum  speed:  ~2  m/s.  Run  description:  l 
Parallel  uphill,  varying  load  every  10  seconds.  Snow  ? 
condition:  Dry,  hard-packed  natural  snoiv. 


Time  \  s ) 


-  7 


~  -  8 


o  -  9 

Figure  Cl  1 .  File  name:  Field  B4A.dat.  Ambient  tern-  * 
perature:  -10°C.  Maximum  speed:  ~10  m/s.  Run  l 
description:  Parallel  down  short,  gentle  incline  with 
load  of  =N/2.  Snow  condition:  Dry,  hard-packed 
natural  snow. 

-  i; 

o  4  0  80  120  160 

Time  is ) 
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Figure  Cl 2.  File  name:  Field  C2A.dat. 
Ambient  temperature:  -2.0°C.  Maximum 
speed:  =9  m/s.  Run  description:  Parallel  down 
short,  gentle  incline  with  load  of  =N/2.  Snow 
con  on:  Soft,  fresh  snoiv. 
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Figure  C13.  File  name:  Field  B4*.dat.  Ambient  tem¬ 
perature:  -9.0°C.  Maximum  speed:  ~10  m/s.  Run 
description:  Parallel  down  short,  gentle  incline  with 
load  of~N/2.  Snow  condition:  Dry,  hard-packed  natural 
snoiu. 


Figure  C14.  File  name:  Field  B10*.dat.  Ambient  tem¬ 
perature:-!  1 .0°C.  Maximum  speed:  ~10m/s.  Run  des¬ 
cription:  Parallel  down  short, gentle  incline  with  load  of 
-N/2.  Snow  condition:  Dry,  hard-packed  natural  snow. 


J  Figure  C15.  File  name:  Field  C4*.dat.  Ambient 
_  temperature:  -1.4°C.  Maximum  speed:  ~ 9  m/s. 
_i  Run  description:  Parallel  down  short,  gentle  in- 
J  dine  with  load  of  =N/2.  Snow  condition:  Soft, 
-  fresh  snow. 

J 

ICO 


Figure  Cl  6.  File  name:  Field  C10*.dat.  Ambient  tem¬ 
perature:  -1 .5°C.  Maximum  speed:  ~9  m/s.  Run  des¬ 
cription:  Pnridlcl  down  short, gentle  incline  with  load 
of  =N/2.  Snow  condition:  Soft,  fresh  snow. 


Figure  Cl  7.  File  name:  FiehlAA.dat.  Ambient  tem- 
pernture:-9.8°C.  Maximum  speed:  =6  m/s.  Run  de¬ 
scription:  Parallel  down  short,  Aentle  incline  with 
load  of  =N.  Snow  condition:  Pi  y,  hard-packed  natural 
snow. 


Figure  CIS.  File  name:  Field  A35.dat.  Ambient  tem¬ 
perature:  -6.7CC.  Maximum  speed:  =6  m/s.  Run  des¬ 
cription:  Parallel  down  short ,  gen  tie  incline  with  load  of 
=N.  Snow  condition:  Dry,  hard-parked  natural  snow. 


Figure  C29.  File  name:  Field  813.dat.  Ambient  tem¬ 
perature:  -2>yC.  Maximum  speed:  ~6  in/s.  Run  de¬ 
scription:  Three  repetitions  of  paralleling  down  one- 
third  of  short,  gentle  incline  with  load  o/=N.  Snow 
condition:  Dry,  hard-packed  natural  snow. 
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Figure  CIO.  File  mime:  Field  B4C.dat.  Ambient 
temperature:  -1 1 .6  'XT.  Maximum  a  peed:  =  70  m/s. 
Run  description:  Parallel  down  short,  gentle  in¬ 
cline  with  load  of  =N/2.  Snow  condition:  Dig,  hard- 
packed  natural  snow. 


100 


Figure  C21.  File  name:  Field  41U.dat.  Ambient  tem¬ 
perature:  -9 .5  CC.  Maximum  speed:  -6  m/s.  Run  des¬ 
cription:  Parallel  down  short,  gen  tie  incline  with  load 
of  =N.  Snow  condition:  Dry,  packed  natural  snow. 


A 

~i 

_ !  FigurcCII.  Filename:  Field  A17.dat.  Ambient 
temperature: -9. 3CC.  Maximum  speed:  m/s. 

—  Run  description:  Parallel  down  short,  gentle 
-  incline  with  load  of  -N,  followed  by  varying 

—  load  up  Poina  lift  with  constant  speed  of  ~  2  ml 
-j  s.  Snow  condition:  Dry,  hard-packed  natural 


Figure  C23.  File  name:  Field  4A.dat.  Ambient  tem¬ 
perature:  -9.5  CC.  Maximum  speed:  -  6  m/s.  Run  de¬ 
scription:  Parallel  down  short,  gentle  incline  with 
load  of  - N/2 .  Snow  condition:  Dry,  packed  natural 
snow. 
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Figure  C24.  File  inline:  Field  4A.dat  (nnsiiwotlied). 
Ambient  temperature:  -9.5  °C.  Maximum  speed: 
m/s-  Run  description ■  Parallel  down  short,  gentle 
incline  with  load  of  ■= N/2 .  Snow  condition:  Dry, 
packed  natural  snow. 


APPENDIX  D:  PLOTS  GENERATED  BY  THE 
FINITE  ELEMENT  COMPUTER  MODEL 

The  transient  response  from  5  seconds  to  steady  state  for  case  1 .  Parameters:  Ambient  temperature  T0 
=  -5°C;  convection  coefficient:  75  W/nv  0°C.  The  magnitude  and  location  of  the  maximum  and  n  mimum 
temperatures  are  indicated  to  the  right  of  the  plot. 
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Figure  Dl.  Transient  response  at  5  seconds. 
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Figure  D2.  Transient  rcsfmiscat  10 seconds. 
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Figure  D3.  Transient  response  at  17  seconds. 
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Figure  D4.  Transient  response  at  17  seconds. 
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Figure  D5.  Transient  response  at  57  seconds. 


Figure  D6.  Transient  response  at  127  seconds. 
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Figure  D7.  Transient  response  at  227  seconds. 
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Figure  D8.  Transient  response  at  377  seconds 
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A  facsimile  catalog  card  in  Library  of  Congress  MARC  format  is  reproduced 
below. 
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